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ADDENDUM 


The results of Volume 2 were discussed recently with Professor 
John Breakwell of Stanford University. His comments regarding 
the matching (Section A14) pointed up a deficiency in the analysis. 

In order for the results of Section A17 to be valid to order p^.the 
guage function discussed in Section A14 should be p rather than 
This stronger guage function would require the inclusion of 
higher order singular terms in both the outer and inner expansions 
in order to satisfy the limit defined by Equation (A14-1). The 
inclusion of such singular terms is somewhat laborious and has 
no effect on the results obtained from the matching and the 
subsequent initial and boundary value solutions derived from these 
results. The matching of such terms will be discussed in a 
forthcoming paper by Breakwell and Perko. 

The limits on the overlap domain depend on which singular terms 
are included in the matching. The overlap domain 2/5< a < 1/2 
found in Section A14 is valid but not unique as other combinations, 
such as 1/2 < a <3/5, are possible. An overlap domain which 
includes a = 1/2 is also possible but requires the inclusion of more 
singular terms than does either of the other overlap domains 
mentioned. The actual choice of overlap domain also has no effect 
on the results of the matching. 


ll 
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INTRODUCTION 


This report contains the analytical derivations of the second order asymptotic 
boundary value solutions for lunar and interplanetary trajectories which have 
been formulated under Contract No. NAS9- 10526 for the NASA Manned Space- 
craft Center. It is a supplementary document to the final report of the 
contract study. * Whereas the final report presents only the results of the 
study effort, this document contains the step by step derivations of both the 
initial value solution to the problem of N-bodies and the boundary value solu- 
tions which were designated in the contract work statement. 

The analysis is divided into three sections. Section A contains the derivation 
of the second order asymptotic solution starting from the differential equations 
of motion for the N-body problem. It includes derivations of both the outer 
and inner expansions, their behavior in the overlap domain, a detailed dis- 
cussion of the matching, and the solution to the initial value problem. 

Section B contains the derivations of a number of boundary value problems 
for both lunar and interplanetary applications. These derivations are based 
on the results of Section A. 

Finally, Section C contains discussions of how the two main types of perturba- 
tion terms in the asymptotic solution are evaluated numerically. This 
includes explicit formulas for evaluating the well known linear state transition 
matrix. 

The derivations are given primarily by showing the mathematical steps 
involved. A minimum amount of discussion is presented in this document. 

For expanded discussions (and a corresponding minimum number of mathe- 
matical expressions), as well as numerical results showing the accuracy and 
computation speed of the asymptotic solutions, the reader is referred to 
Reference 1. 


I 


The notation used in this report is a combination of that of Lancaster and 
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Carlson. In general, each parameter is defined as it is introduced but some 
which have only mathematical meaning and serve an intermediary role are 
defined only by an equation. Scalers are written as x or X and vectors as 
x or X. A matrix G(x) and a tensor H(x) are also used. In addition, a bar 
over a parameter indicates that it applies specifically to an inner solution. 
Finally, the order of a particular term in an expansion is given by the 
exponent of the parameter p which precedes the term, i. e. , p. n is order n 
or O(n). 



Section A 


SECOND ORDER ASYMPTOTIC SOLUTION TO THE PROBLEM OF 

N- BODIES 


A 1 N-BODY EQUATIONS OF MOTION 

The problem of N-bodies will be defined as follows: The motion of a small 
body of negligible mass is to be determined subject to the gravitational 
forces of a primary body of mass m Q and N-2 secondary bodies of mass ’ . 
nx, i = 1, 2, . . . , N-2, where nx « m^ for each i and the motion of the N-2 
secondary bodies relative to the primary body is assumed to be a known 
function of time. Denoting the position of the small body relative to the 
primary body by r* and the positions of the secondary bodies by pL : the 
differential equation for the small body motion is 




-G 


N-2 



m. 

l 



(A 1 - 1 ) 


where G is the gravitational constant. 
Now define dimensionless variables 


r = r /L 
t = t / T 


whe re 


L = a^ = semi-major axis of the motion of the j body 


th. 

2uT = P. = period of the j body motion 


. (A 1-2) 
(A 1-3) 


(A 1-4) 
(A 1 - 5) 
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Also define the dimensionless mass ratios 


: : i 


M- 


. = m. / m 


1 o 


(Al-6) 


Then (Al-1) becomes 


d 2 r r V 

~7T £ v-i 

dt r i=l 


N-2 


r - p. p. 

~ + _ZL 

3 3 


£- £i 


p i 


(A 1-7) 


(Al-7) is the dimensionless differential equation of the small body where the 
unit of length is now L* and the unit of time T * . Defining 


f(x) 


- x/x 


X 


d 2 x 

dt 2 


(A 1- 8) 
(A 1-9) 


reduces (Al-7) to 


r = f_(r) + F(r, p.) 


(A 1- 10) 


where 


. N-2 

ESL’ £i> = 2 p i " -Ei* + - ( -Ei } ] 


(A 1 - 1 1 


In (Al-10) r_ and f_(£^) are both of order unity (zero order) while .F is order p. 

(first order), except in the exceptional case where r - p. = 0 (p. ) when F is 
_ 1 — —l l 

order p^ . However, in the analysis of transfer trajectories the exceptional 

case is the case of interest since transfer trajectories have at least one-close 

approach to a secondary body. The change in the order of F indicates that 

4 

(Al-10) represents a type of singular perturbation problem . This assump- 
tion will be verified by developing a solution to (Al-10) by the method of 
matched asymptotic expansions. 
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In addition to (Al-10) the equations of motion of the secondary bodies will be 

th 

of interest. For the k body they are 

N-2 

£ k = (l+n k )l<p k )+ 5! ^ [i(£ k - Ei) +l(P i )l (Al-12) 

i= 1 1 J 

i^k 

Al. 1 Interplanetary Equations of Motion 

For interplanetary trajectories the primary body is the sun and the secondary 

til 

bodies are whatever planets one wishes to include in the model. The j 

planet used to determine the length and time scales, L* and T*, will normally 

be the target planet of the trajectory. All of the p. 's will be small, the 

-3 1 

largest being about 10 for the planet Jupiter. In this case it is obvious that 
IT is small in (Al-10) except when the trajectory makes a close approach to 
one of the planets. 

• 

In this case the coordinate system is centered at the sun and it is the motion 
of the sun about the solar system center of mass which produces the f_(jD. ) 
term in (A 1- 1 1 ). 

A1.2 Cislunar Equations of Motion 

In cislunar space the primary body is the earth and the secondary bodies are 
the moon and the sun. The j^ body will be the moon so that L* is the semi- 
major axis of the moon's motion about the earth and 2-rrT* is the moon's 
period. The coordinate system is centered at the earth and the fjjx) term in 
(A 1—11) is due to the motion of the earth about the center of mass. 

(Al-10) becomes 

r = £(r) + p M [f_(r - p M ) + f_(p^)] + p g [f_(r - ^ ) + fj£ g )] (A 1-13) 

Since 

J , _ , 

= m S /m E 
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it is obviously not small. Let 


V- = M- 


M 


M s = ^s 7 ^ 


(Al-14) 
(A 1- 15) 


and (Al-12) becomes 


r = f (r ) + p 


[i<£ - P M ) + M s[l ( r - P S } + I ( £s ) J 


(A 1- 16 ) 


It is necessary to show that the term proportional to is actually order p, 
i. e., to show that the sun's effect is the same magnitude as the moon's. Let 


> = a p 

-s s — s 


(A 1- 17) 


where a is the mean distance of the sun from the earth (in dimensionless 
s . 

units). Then 


r - p = -a (p - r/a ) 
— J-s s — s — , s/, 


(A 1- 18) 


and, since r«a 


IS* ~ £ a > = " iS£ a " £ /a s ) /a s 


4 i'£s> + 4 C(_P s )r + 0(4-) 
a a a 

s s s 


(Al-19) 


where G(£ ) is the gravity gradient'matrix defined in Section A 2. Also 
s 


i<p s ) = i(P s )/a^ 


(A 1-20) 


giving 


f(£ - p g ) + f(p s ) 


~J G(p g )r +0(-^) 

a a 

s s 


(A 1-21) 
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or 


M 


[f(£-E s )+t(£ s )] 


= o 


^ M s G (P s )r\ 


(A 1-22) 


In (A 1-22) 

M s = M- s / m- = ( m s/ m E)/* m M /m E) = m s /m M = 2. 7 x 10 ? (Al-23) 

G(P ) = 0(1) (since p = 0(1)) 1 (Al-24) 

— s — s 

r = 0(1) (A 1-25) 

a^ = a /L ^ = (lau/0.0026 au)^ = 5. 7x10^ (A 1-26) 

s s 

Therefore 

M s [i(£ " £ s ) +f(£ s )] = (M g /a*)xO(l) = 0.47x0(1) 

= 0(1) (A 1- 27) 

and 

M- m s [i(£ - £ s ) +i(£ s )] = 0(n) (A 1-28) 


In later analysis (Al-17) will not be introduced and it will be implied that 
(Al-27 and 28) hold for cislunar space. Thus the sun and moon are expected 
to have comparable effects on cislunar trajectories. 


A2 EXPANSION OF f_(x) 
Let 


f_(x) = f.(x) e. 


(A2 - 1) 
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where f are components of f and e^ are orthogonal base vectors. Also 
let 


x = x + 6x 


(A2- 2 ) 


where 


2 3 

6x = ex j + e + 0(e ) 


= (exj. + £ x 2i + . . . ) e. 


(A2-3) 

(A2-4) 


Substituting (A2-2) into (A2-1) and expanding in a Taylor series gives 

9f. (x ) . 9 2 f . (x ) , 

f. (x) = f . (x ) + — — 6x. + ■=■ —k — * — 6x. 6x. + 0(6x ) 

l— l—o 9x. i 2 9x.9x, i k ' ' 

J J J k J 


9 f .(x ) r 

f i ( ^o )+ “9^— L“lj 

J . 


+ e Z x 2j + 0(e 3 ) 


. 9 f. (x ) - , 

, 1 i o' 2 , 3 X 

"z ix.av e x ij x ik + ° (£ > 

J k 


(A2 - 5) 


Let 


G..(x) = 

ij - 


9f.(x) 

^T 


9 

9x. 

J 



x. 


jlLA 

i 9x. 1 3 J 3 9x. 

J\ x / x J 


(A2- 6) 


But 


x 


x x 
m m 


. 2x 3x = 2x 


9x. 
J 


9x 


m 


m 9x. 
J 
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(A2- 7) 


9x 

9x. 

J 


x 9x 
m m 


x 


9x. 

J 


and 


9x. 

0 — = 6 .. = kronecker delta (A2-8) 

Xj lj 


Therefore 


9x _ X j 

9x. x 
J 


Using (A2-7) and (A2-9) in (A2-6) gives 


(A2-9) 


3x.x. 6 .. 

J XX 


(A2^10) 


Now let 


Sijkfe) 


3 2 £.(x) 

9x.9x, 

J k 


x.x. 6 . 


a G..(x) 
9x k 




15x.x.x, , 

-Aji— + - 5 (x. 6 jk + Xj 6 . k +x k 6 i . ) 


(A2- 11) 


Summarizing, if 


f(x) = -x/ 


x 


(A 2- 12) 


and 


x 


X + £ X 

— O — 


+ £^ x^ + O(e^) 


(A2-13) 
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then 


f_(x) = MXq) + sG^) xj + e 2 [G(x o )x 2 +^H (^xj] +0(e 3 ) (A2-14) 


with 


G..(x) 


3x.x. 

1 J *J 

5 3 

x x 


(A2-15) 


H. (x) = 


15x.x.x, ^ 

^ + ^T ( x ; 6 ^ + x ; 6 ;u + x iAJ 


—ijk'—' v 7 ' ^5 jk j ik ~k”ij 


A2-16) 


x x 


The matrix G(x) is the gravity gradient matrix and is also defined by 


G(x) = 


_ df(x) 


dx 


(A2-17) 


An expansion of G(x) is similar to (A2-5): 

9G. . 9 2 G..(x ) - 

G. .(x) = G. . (x ) + a ^ (x ) 6x, + -3 — U — — 6x, 6x + 0(6x ) 
ij — ij '—o' 9x^ v — o k 3x^9 x^ k l 


3G..(x ) ? 

= G ij ( S 0 )t bi ° [“ij + E x 2j + ° <E >J 


, arc.. he. ) , 

1 n'— o 2 




(A 2 - 1 8 ) 


Let 


T 

~ ijkf 


9 G..(x) 

LLzl 

3 V X « 


= -J- H.., (x) = 

dx „— iJ k ~ 


105 x.x.x, x. 
i ] k I 

3 


15 

~T 

x 


h (x j 


6 , . - x„d., ) 
kf £ jk 


+ x.(x. 6 . . - x. 6 .. ) + x, (x. 6 . - x 6 ..) 

j' k if i ik' k i j£ £ ij _ 


x 5 ( 6 i? 6 jk + 6 jf 6 ik + 6 k? 6 ij } 


(A 2 - 19) 
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Then 


G(x) = G^) + eH(x Q )x 1 +£ 2 [h(x q )x 2 + T(x q )Xj/ 2] +0(e 3 ) (A2-20) 

A3 LINEAR DIFFERENTIAL EQUATIONS; STATE TRANSITION MATRIX 
In order to develop an analytical solution to (Al-10) it is necessary to solve 
a system of first order differential equations of the form 

(A3 - 1 ) 
(A 3 - 2 ) 

where G(t) is the gravity gradient matrix and the F (t) are given functions of 
time. This system can be written in the condensed form 


r 

— n 


v 

— n 


-n 


G(t)r + F (t) 
— n — n 


x =. K(t)x + u (t) 
— n- ' — n — n - 


(A3-3) 


whe re 


x 

— n 


-n 


(A3-4) 


u 

— n 


O 

F 

. — n. 


(A3 - 5) 


K = 


<0 I 

iG O) 


(A 3 - 6 ) 


3 5 6 

The solution of (A3-3) is well documented ’ ’ and can be written in the form 


x (t) = 4>(t, t ) x (t ) + 

— n o — n o 


t 

j $(t,T) U n 


(T)dT 


(A3- 7) 
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where the matrix $ satisfies the differential equation 


d$ , > 

dF (t > T) 


K(t)$(t,T), $(t, t) = I 


or 


d<$ . . 

37 (, - T) 


-$(t,T)K(T), $(T,T) = 1 


(A3- 8) 


(A3- 9) 


It is easily verified that (A3-7) is a solution of (A3-3) by differentiation: 


- 10 ) 
t 

o 


X (t) 

— n 


* {t ’ ‘o n^o ' 


f $(t, T )u n 


(T)dT + $(t, t)U n (t) 


(A3 


The last two terms in (A3- 10) come from differentiating the integral in 
(A3-7). Replacing $(t, t Q ) and <$(t, t) by (A3-8) gives 


x (t) 
— n 


K(t)$(t,r)x n (t o ) + 



K(t)$(t,T)U n (T)dT + U n (t) 


= K(t) 


${t > T) in (t o 


) + 



4 , (t. T )U; n ( T ) dT 


+ H n (t) 


(A3- 11) 


Replacing the bracketed term in (A3- 11) by (A 3-7) gives 

x n (t) - K(t)x n (t) + u n (t) (A3- 12) 

and (A3-12) is identical to (A3-3). QED. 

3 5 

The matrix $ is the well known state transition matrix ’ , also known as the 

6 7 

matrizant and the fundamental matrix . It is a 6 x 6 matrix which can be 
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partitioned (following Carlson's notation) into four 3x3 partial derivative 
matrice s 


(A3- 13) 


(A3- 14) 


A(t ' T > 

= C(t,T) 

(A3- 15) 


= D(t,T) 

(A3- 16) 


= G(t) A(t,T) 

(A 3 - 1 7) 

a 

07 D(t,T) 

= G(t) B(t,x) 

(A3- 18) 

37 a <*- t > 

= -B(t.T) G(t) 

(A3- 19) 

af B (t,T) 

= -A(t,T) 

(A 3 -20) 

57 c < l ' T > 

= -'D(t.r) G(T) 

(A 3-21) 

37 d "' t > 

3 

o 

= -C(t,T) 

(A3-22) 

A(t,t) = 

A (T, T ) = I 

(A3 - 23 ) 


$(t,T) = 


/3r(t) 

ar(t)\ 

3r(r) 

3v(t) \ 

. 3i(t) 

3v(t) I 


3v(t)/ 

A(t,T) 

B (t, T )\ 

\C(t,T) 

D(t,T )j 


The four partial derivative matrices have derivatives given by" 
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B(t, t) = B(t,t) = O 

(A3-24) 

C (t, t) = C(t, t ) = O 

(A3-25) 

D(t, t) = D(t,t ) = I 

(A3-26) 

From (A3- 15) through (A3-26) various Taylor 

3 

derived. Some examples are 

series expansions can be 

A (t, t) = I + G(t)(t-T) 2 /2j + O(t-r) 3 

(A3-27) 

B(t, t) = I(t-r) + G(t)(t-T) 3 /3: + O(t-T) 4 

(A3-28) 

C(t,T) = G(t)(t-r) + O(t-T) 2 

(A3-29) 

D(t, T ) = I + G(t)(t-T) 2 /2: + O(t-r) 3 

(A3- 30) 


Some additional forms of the partial derivative matrices are given in 
Section C . 

Two properties of the state transition matrix which are used in later 
sections are 

Gftg.tj) = $(t 1 ,t 2 )' 1 (A3-31) 

■■ ®(t 2 .tj) = $(t 2> t o )$(t o> t 1 ) (A3-32) 

These properties are especially useful in formulating the solutions of 
Section B. 

A4 OUTER LIMIT 

The outer limit of (A 1-10) is the limit where r - jx = 0(1) for all 
i = 1, 2, . . . , N-2. In other words, the outer limit specifically excludes 
the exceptional case where £ - = O(p^). Therefore the function F in 
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(A 1-10) is always order in the outer limit. In the outer domain, i.e., 
the domain defined by the outer limit, the solution of (A 1-10) is assumed 
to be of the form 


+ |jl r 


-1 


-2 




-3 


where 


(A4- 1 ) 


\ 4 - 2 ) 

(A4-1) is a representation of the solution r in the form of an asymptotic 
expansion in powers of p. It is an exact solution of (A 1-10) when p = p^ = 
p^ = . . . = p^ 2 = 0 and approximates the exact solution as long as p is small. 


A5 OUTER DIFFERENTIAL EQUATIONS 
Substituting (A4-1) into (A 1-10) gives 


. . . . . . 2 . . 3 

L = Lq + P£i + £ 2 + °(l J - ) 


f_(£) = fjr^) + + p 

N-2 


G (r )r_ + — H(r )r 2 
v —o —2 2 i— '—o—l 


+ 0(p 3 ) 


F(£> £i ) = 2 M- M i [i(£o - £i) + £i 

i _ 1 L 

+ 0(p 2 ) + f (£.)] . , . 


whe re 


(A 5- 1 ) 
(A5-2) 


(A 5-3 ) 


M. = p./p 

i ~i ~ 


(A5-4) 


Equating power of p in (A 5-1) through A 5- 3) gives 


r = f(r ) 
— o — — o 


r. = G(r )r, + F,(r , p. ) 
— 1 '—0—1 —l—O -s-l 7 


(A5- 5) 
(A 5 - 6) 
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(A5-7) 


—2 G ^^o } — 2 + — 2 ( ^o’ 'v P i ) i 

. • ■ ' * * ; 

or, in general, 


— n ' G <Io>£n + ^n ( ic £n-l' £i> 


(A5- 8) 


The functions F ^ and are given by 

•' L ( £o- £i> - j M i [!(£<,- £i> +K£i)] 

i=l L 


(A 5-9) 


EzL,-. £ lf Hi> 


N - 2 

|h(£o>£i + Z M i G L 

1 = 1 


£i>£l 


(A5-10) 


The general term would involve tensors up to order (n + 1), Because of 
this behavior it appears somewhat impractical to go beyond n - Z which 
already includes the third order tensor Hjr^). 

The differential equations (A5-5), (A5-6) and (A5-7) are the zeroth, first and 

■ ; . , . • i 

second order outer differential equations. The equation for is simply the 
two-body differential equation while those for r^ and r^ are of the linear type 
discussed in Section A3. 


A6 OUTER SOLUTION 

The zeroth order outer differential equation can.be written 

,2 

dr r 

-o — o 



(A6-1) 


The solution to this differential equation is well known and many forms exist. 
The motion is elliptical with respect to the primary body and a useful form is 


Io (t) 


£ oV 


r (t ) 
— o o 


g D (t) 


(t ) 
> o 


(A 6- 2) 
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where r 


is the two-body position and v the velocity 


dr 


— o 



(A 6- 3) 


The functions f and g are infinite series in time but have closed form 
o "o 

expressions as functions of eccentric anomaly E where 


n (t - t ) = E - e sin E 
o po o 


(A 6- 4) 


In (A6-4) n Q is the two-body mean motion, e Q the eccentricity and t the 
time of pericenter passage. The functions f and g are given by 


f (t) = 1 - a 1 - cos AE (t 

o , o L 


,t )~|/r (t ) 

o J/ o o 


(A6- 5) 


g (t) = ("r (t ) sin AE (t, t ) /( n a ) 

"o' o' o' o'J/ ' o o 


+ e sin E (t ) 1 - cos AE {t, t ) / n 

O o L o J/ o 


(A6-6) 


where 


AE(t,t o ) = E(t) - E(t o ) 


(A6-7) 


The higher order solutions are obtained by first writing (A5-8) as 


r = v 
— n — n 


(A6-8) 


v = G ( r ) r + F (r , . 
— n '—o' — n — n — o 


•• in-1' £i> 


(A6-9) 


Since (A6-8) and (A6-9) are similar to (A3-1) and (A3-2) the solution follows 
from (A3-7) 



®<t, t Q ) 


r (t ) 
— n o 


v (t )j 
l— n v o • 


c 




0 




(A6-10) 
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Using (A3- 14) gives 


£ x (t) 


A(t, t Q ) £ 1 (t Q ) + B(t, t Q ) + 


J B(t,r) FjfrJdT 

t 

o 


(A6-11) 




A(t, t Q ) 






£2 


(t ) + 


/ 


B(t, t ) F 2 ( T )dr 


(A6-12) 


The outer solutions are therefore given by (A6-2), (A6-11) and (A6-12). They 
are functions of the time t, the initial time t , and the initial position and 
velocity 




^o ( V + ^l (t O ) £-2 (t Q } 


V (t ) + [iV,(t ) + U v_ (t ) 
— o' o 1' o' r —2' o' 


(A6-13) 

(A6-14) 


It will be shown that using these solutions the outer expansion (A4-1) contains 
a non-uniformity when _r - jx = 0(|ju). That is, the individual terms in the 
expansion do not remain small compared to the preceding terms. In this 
case it is necessary to investigate another limit of (A 1-10). 


A 7 INNER LIMIT 

When the trajectory representing a solution of (Al-10) passes close to one 
of the secondary bodies the outer solution is no longer valid due to a non- 
uniformity in the outer expansion (A4-1). In order to study another limit the 
origin is transferred to the kth secondary body and the length scaled such that 
the new position vector is order unity. Such a transformation is given by 


—a 


<£ - Ek^k 


(A 7-1) 


It is obvious that when x_ - = 0(p.^) then x_ Q = 0 (1). It is also necessary 

to make a similar transformation of the independent variable, i. e. , 


“ (t " t pk ) ^k 


(A 7- 2) 
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where t ^ is some fixed time, such as time of closest approach, associated 

with motion close to p, . When f 

— k 


;n t - t k = 0(|i^) then tp = 0(1), 


The next step in determining the inner limit is to transform (A 1-10) into the 
new variables. From (A7-1) 

„ . ■ ... i . , • i' 


£ = £k + -a 
and differentiating gives 


• • 

£ = £k 

• • • • 

£ = £ k 

Also 

i<£> = HEfc + (*£££,> 

and, using (A2-14) 


dr 

+ 

^k dtp 


+ M- 


, 1 , d 2 r 

a- 2p —a 

dt I 



I 7 rn 7 'Kn 

L (£) = i(£ k ) + c (£ k ) £n + P k H(£ k )r a /2 + °^k ) 


(A 7- 3) 


(A7-4) 


(A7- 5) 


(A7- 6) 


Likewise 


£(£.£i> 


N-2 

“iP^i’-k 1 £a - Ei> ‘IL’j 

1= 1 L -> 


r t ¥ r 

*k IK£«> + f( £k> + 4/i 1 

J 1= 1 L 
i^k 

+ ^G(£k-£i>il O + 0 <¥> 


(p k -■£•)+ l(£i) 


(A 7- 7) 
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where 


LK ^ 


HL a )/v- k 


(A 7- 8) 


Substituting (Al-12) into (A7-5) and (A7-8) into (A7-7) and then substituting 
(A7- 5), (A7- 6) and (A7-7) into (Al-10) gives 

N-2 

I(£ k ) + P k K£ k > + E M-i[l(£ k - £j> + K£i)] 

i^k 

l(£ k ) + p£ G ^L a + I P k * *£(£ k } £« 

+ n k " 2a £(£ k ) + P k £(£ k ) + °(P k & ) 

N-2 

+ E Pi [l(£ k - £ i )£(£ i ) 

1= 1 L 

i^k 

2a ~\ 

+ P k G (£ k ' £i> £a + 0(P k )J (A7-9) 


+ P 


g-2p —a 
2 

ndto ■ 


or, cancelling similar terms gives 


d £« 

dtp 


= P 


l+2p-3o 


,2p 


1 a+ 2p 


*<£«> + P k ^ C( £k )r 0 + | ^ H( £k ) £ ; 


N-2 

+ P k P E P i G( £k - £i )r ff + 0( h x k "^, nj^p) 

i= 1 


2U+2P u+2p 


(A 7- 10) 



(A7-10) still includes functions of t through^ and jx. From (A7-2) 


t 



^k fc f3 


(A7-11) 


Let 




(A7-12) 


giving 

t = fc k + (t p + T k^ (A7-13) 

The significance of t^ and will be demonstrated later in the matching. 
Using (A7- 13 ) 


G (E k (t)) = G (£k (t k } ) + ^k ^ G ( £ k {t k ) ) (t P + T k ) + °(^k P ) (A7-14) 


The derivative dG/dt is found as follows: 
F rom (A2- 15) 


dt G ij^ 


d 

dx, 


/ 3x.x. 

6 - \ 

dx, 

1 5x.x. 

9x 

I 1 1 

U 

k - 

1 J 

V 5 ' 

V x 


1 dt 

6 

X 

9x. 

k 




Using (A2-9), (A7-15) becomes 


3x. 36.. a ' 

_j , _i] dx 

+ 5 5 ik 4 9x. 

x x k 


x, 


(A7-15) 


ft G ij<2> 


15x.x.x, , 

+ — =• (x.6-, + x. 6., + x, 6..) 

n 5 i jk j lk k lj 


x 


x 


Sijk'-’-k 


(A7-16) 
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Then 


C(p k (t)) = G k + p£ H k p k (t k ) (t p + r k ) + G(pf) 


(A7-17) 


where 


G k = G (£k (t k>) 

H k = S(p k <‘k>) 

Putting 

G (E k (t k ) ' Ei (l k>) = G k 


(A7- 18) 
(A7-19) 


(A7-20) 


(A7- 10) can be written 


d 2 r 


dt2 = 11 k + ^ - ( -o } + ^k G k -a + ^k — k -Ek (t k ) -a (t p + T k } 


N- 2 

. 1 a+ 2(3 „ 2 , 2PST _i 

+ 2^k P H k r^+.p k ^ ^ G k r &; 

i=l 

i^k 


+ o(p 2 ^ 


(A7-21) 


The inner limit is partially defined by balancing the inertial term with the 

t h. 2 2 

gravitational term of the k body, i, e. , matching (d £ a /dt ) with f_(r ff ). The 
two terms will be balanced if both are order unity. This is obtained if 

1 + 2(3 - 3a = O (A7-22) 
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Also, from (A7-4) 


dr 

O' 

dtp 



(A 7-23) 


In (A7-1) it is assumed that the numerator, r - is small; this gives a 
close approach to the kth body. In (A7-22) it is not necessary for _r - p^, the 
relative' velocity, to be small. In fact, this velocity difference will depend on 
the energy of the outer solution and the only conclusion that can be made at 
this point is that it is order unity. Thus it is most reasonable to make the 
denominator in (A7-23) the same order, i.e., let 

a - p = O (A 7- 24) 

The inner limit will be defined by (A7-22) and (A7-23), i.e., 

o = P = 1 (A 7- 2 5 A) 


As a result (A7-21) becomes 


d 2 r 


~a 


dt 2 p 




— k£ k *V -a (t p + T k ) 


N--2 

+ H, rl + M.G? r 
2 — k— a *—> i k —a 

i= 1 

i^k 


+ O(JJL^) 


(A7-25B) 


In the terminology of singular perturbation theory, (A7-24) is called a 
distinguished limit. 


The inner variables will be defined by 

R k = (£-.£ k >/n k 

S k = (t - t pk )/ ^k 

24 


(A7-26) 
(A 7- 27) 



and the associated differential equation is 


d 2 R. 

* = i(R k ) + P(R k , £ k . £i ) 

dS, 


(A 7- 2 8) 


whe re 


P(R k - £k» £i } 


G k ^k + ^k 


-k £k ( V -k (S k + T k^ 


N-2 


+ 4 H. R? + A M. G* R, 
2 — k — k i k — ; 


i= 1 

i^k 


+ 0(p*) 


(A7-29) 


3 

The term of order p^ in P is included here although it is not actually needed 

in a general second order theory. It is important in problems where only an 

inner solution is used. It may also be important in problems where the 

model includes satellites of secondary bodies, i.e., including moons in the 

th 

analysis of interplanetary trajectories. Suppose the m body is close to the 

th i ' *•••;:. ’ / 

k body such that 




•k — m — m 


£ P. 


(A 7-30) 


Then 


G IT = C <£k-£m> 


= G(ep ) 
' — m 


G(p ) / e* 
— m 


(A7-31) 


Then the te rm 


3 A . ^m 

i M G i 

k m k 


(p? p /e 3 ) G { p ) 
r k ^m — m 
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For the moon, earth, sun case 


a = m /m = 3.7x10 
m o 


-8 


so that 


_ 3 

e = 2. 6 x 10 au 


<“m /£ > 


2.1 = 0 ( 1 ) 


and 


[±l M 


G™ 
m k 


O(M-u) 


(A7-32) shows that the effect of a nearby moon e. , |i 


? M g“ R 
k m k 




(A7-32) 


on the 


inner solution of an interplanetary trajectory may be the same order of 

2 

magnitude as the effect of the Sun (i.e. , R^)* 


In the inner domain, i.e., the domain defined by the inner limit, the 
solution of (A7-28) is assumed to be of the form 


R k 


= *ko + 


M-k R 


■k2 


+ ^k — k3 + ^k — k4 + 


(A 7-33) 


(A7-33) is also an asymptotic representation and reduces to the exact 
solution as |i^-^0. 

It should be noted that another distinguished limit exists besides (A7-24). If 
the velocity difference in (A7-22) is assumed to be smaller than order unity 
then the other limit is 


O' = 1/3, (3=0 


in which case (A7- 10) becomes 


d ^ r 

-f* = Kl )+G(£ k )r c +0(n 1/3 ) 
dtp 


(A7-34) 


(A 7 - 3 5 ) 


26 



(A7-35) is a form of Hill's equation and it's solution is beyond the scope of 
this study. 


A8 INNER DIFFERENTIAL EQUATIONS 
Substituting (A7-33) into (A7-28) gives 


d 2 R 
— k 

dS? 


d -ko , 2 d — k2 , 3 d — k3 , ^, 4, 

dsf K dsf K dsf k 

k k k 


K5k> = '(Sko> + ^k °®ko> *k2 


+ >4 G ' — ko ' — k3 + °<4> 


E'-k’ Ek' Ei> = 4 G k 1 ko + 4 


H k i k (t k )R ko (S k + r k ) 


N-2 


♦ki* .2, M i G l 


i=i 

i^k 


+ O(n^) 


Equating powers of in (A8-1) through (A8-3) gives 


d 2 R 


ko 


dS, 




d k k2 

-^r = G( ^ko> « k2 + ^2 <5ko> 


d k k3 

dS 2 


G(R ko ) R k3 + p 3 (R ko ) 


(A8- 1) 


(A8-2) 


(A 8- 3) 


(A8-4) 


(A8- 5) 


(A 8 - 6) 
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where 


-2 ( -ko ) G k-ko 


(A8- 7) 


L'^ko' = ? «k*ko + «kik'‘k»5:kc < S k + T k» 

N-2 

+ S M i G k R ko <A8 - 8 > 

i= 1 

i^k 

The differential equations (A8-4) and (A8-6) are the zeroth, second and third 
order inner differential equations. The equation for R. is simply the two- 
body differential equation while those for R^ anc ^ — k3 are ^ near type 

discussed in Section A3. 


A 9 INNER SOLUTION 

The solution to the zeroth order differential equation is hyperbolic motion 
with respect to the k secondary body and can be represented by 

— ko < S k> = 7 o ,S k> — ko < S ko> + L < S k> ^ko < S ko> (A9 ‘ 


The functions f Q and g Q have closed form expressions as functions of the 
eccentric anomaly where 

"k s k = ? k sinh F k - F k (A9 - 2) 


In (A9-2) n^ is the two-body mean motion and e^. is the eccentricity. The 

functions f and g are given by 
o o 


f 

o 


- a k[ l - c ° shAF k< S k' S ko>]/ R ko< S ko> 


(A9-3) 


s 0 = s k - S ko -[ sinhAF k ,S k' S ko ) - AF k (S k- S ko ) ]/“k 


(A9-4) 
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where 


aF k< S k- S ko> = F k< S k> - F k ,S ko> 


(A9-5) 


It was shown by Carlson that it is desirable to choose the initial conditions 
for the higher order solutions to be zero at pericenter of , i.e., the 
higher order solutions vanish at S^. = 0. A slightly more general approach 
is to put 


?:k2 (S ko ) = X-(S,J = 


— k2 v ko 


— k3 ko 


— k3 ko 


(A9-6) 


so that the higher order solutions vanish at . Since (A8-5) and 

(A8-6) are identical in form to (A5-6) and (A5-7) the solutions will be 
similar to (A6-11) and (A6-12) 


— k2 (S k } = f k B(S k* (r) —2 (<r) d(r 

• 4 , 

ko 


(A9-7) 


/' 


— k3 (S k ) = / B(S k* <r) —3 (tr) dcr 

“S. 

ko 


(A9-8) 


where B is a partial derivative matrix evaluated along the hyperbola R^. 


The inner solutions are given by (A9-1), (A9-7) and (A9-8). It will be shown 
that using these solutions the inner expansion (A7-33) contains a non- 
uniformity as 


A10 OVERLAP DOMAIN 

The outer expansion has already been given as 

£(t) = r^t) + M.£j(t) + p 2 £ 2 (t) + 0(p 3 ) (A 10- 1 ) 
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Rewriting (A7-26) using (A7-33) gives the outer form of the inner expansion 
as 

i 

r(t) = ]^<t> + R ko < S k > + >4 R k 2 (S k > 

+ R k3 (S k ) + o(njj) (A 10-2) 

These two solutions are valid in the outer and inner domains, repectively. 

b , ... - ... . 

It is assumed that both solutions are valid in an intermediate domain called 

the overlap domain which connects the outer and the inner domains around 
th 

the k secondary body. It will be shown later in the matching that one of 
the conditions for the two solutions overlapping is that for some t = t^ 


io (t k> = £k (t k> 


(A 10-3) 


where t. was first introduced in (A7-12). The time t, is the nominal time 

K K 

at which a close approach to the k body occurs. The actual closest approach 
time is given by (A7-12). 


In (A6-13) r^O^) is the initial position for the two body ellipse ^(t). If t^ is 
specified and an ephemeris is used to determine £^(t^.) then (A 10-3) gives 


the final position r^t^) for the two-body ellipse. The initial and final posi- 
tions and the time interval (t, - t ) constitute the standard Lambert problem, 

K O 

Solution of the Lambert problem completely determines r (t). 


The overlap domain will be defined as the domain where t — ►t^ in the outer 

solution and where |S^ | — - <x> in the inner solution. Thus the behavior of the 

outer solution must be determined in a near neighborhood of the k^ secondary 

body and the bejiavior of the inner solution must be determined in a region 

th 

sufficiently removed from the k body. 
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All BEHAVIOR OF THE OUTER SOLUTION IN THE OVERLAP DOMAIN 


All.l Zeroth Order 

The overlap domain is defined as the domain where t -* t^. It is therefore 
necessary to expand the outer solution in terms of t - t^. The zeroth order 
solution is easily expanded since a two-body ellipse is defined by functions 
which can be expanded in Taylor series' about t = t^. The result is 

Io (t) = ^ok 4 —ok (t " V + ^ok (t ‘ V + ^ok (t ‘ t k ) 

+ e ok (t “ ^ 4 °( (t " V 5 ) (All-1) 

where 


—ok 


b u 
—ok 


^ok 


^ok 


r (t, ) 

—o' k 

L (t k> = ^ (t k> 

L (t k )/2 = tfeo'vV 2 

■r o ( t k )/6 = i(r o (t k ))/6 


= #-(io (t k , )L ( *k > /6 

= c (£o (t k >) ^o (t k )/6 


(All-2) 

(All-3) 

(All-4) 


(All-5) 


^ok = lo ^ 24 

= £ [ G (io<v) ^o<v]/ 24 

= [s^M^k' + ^ioM^ioM]/ 24 ,A1I - 6 > 
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,, OUTER SOLUTION 
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The expressions for c . , d . and e , follow from using (A5-5), (A2-17) 
r — ok —ok — ok 

and (A 7- 1 6). 


Since the motion of the k body is assumed to be a known function of time it 
can also be expanded in a Taylor series: 


E k <‘> = E ko + Ekl" - V + Ek2 (t - l k )2 4 Ek3 (t - ‘k> 3 


4 Ek4 + °(<‘ - V 5 ) 


where 


£ko = £k (t k> 

£ki =. ik (t k> 
EkZ = E-k^ 2 

Ek3 = Ek (t k )/6 


= 5T i(Ek(‘k>) + EEk"’]/ 6 


G (E k ( t k>) £k (t k> 4 ^Ek" 


E k 4 = Ek (t k )/24 


d_ 

dt 


^ifekM ]/ 24 4 °(E> 


“(fikM Ek (t k> 4 G (E k (t k>) t(Ek (t k>) 


24 


(All-7) 


(All-8) 

(All-9) 


(All-10) 


(All-11) 
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+ O(p) 


(All-12) 


and 


/ 

/ 


£k = (X<*> - itfc(‘))] t=tk A 


N-2 

= M ki(£k (t k>) + I - £i (t k>) 

i= 1 

i^k 

+ l(£i< t kl)] (All-13) 


£k” ■= ^[i>k (t) -t(Ek(*k>)] /* 

k 


N-2 


- M k G k + 2 M i _ G k (■£k (t k ) " JMV) 


i= 1 

i^k 


G (£i(t k )) £i(t k )] 


(All-14) 


The expressions for -£k3’ — k4’ £k an< ^ £k ^°^ ow ^ rom using (A 1-12), 

(A2-17), (A5-5), (A7-16), (A7- 18) and (A7-20). 


Now define r. by 
— ko 1 

iko s io - £k 

= ( 2ok - £ko> + ( -ok - £kl> (t - V 

+ <£ok - £k2> (t - *k> 2 + <4>k - £k3> - *k> 3 

+ feok - Ek4» (t - »k> 4 + °( ( * ' V 5 ) (All-15) 
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where 


£ok ‘ -Eko ' io (t k’ ' £k (t k> 

(All-16) 

&ok - £kl = io (t k> - £k (t k' 

(All-17) 

£ok - £ k 2 = [ifcoV) - i(£k (t k>) - ££k]/ 2 

(All-18) 

iok - £k3 = [ c ^o (t k>) £o (t k> - G (£k- (, k>) £k (t k> 


-££k j/ 6 

(All-19) 

Sok - £ k 4 = [5(£o (t k>) 4<V - H (£k ft k>) £k<V 


+ G (w) i(iA>) - G (£k<‘k») L( Ek<v) + 

jzA (All-20) 

From (AIO-3) and (All-16) 


£ck-£ko = ° 

(All-21) 

Let 


^ok " £kl = ^o ( V " £k (t k ) = -k 

(All-22) 

Again using (A10-3) and (All-22) 


-ok " £k2 = "^k^ 2 

(All-23) 

^ok " £k3 = [ G (£k (t k ) ) -k " ^k ]/ 6 

(All-24) 
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£ok - £k4 = [a(£k (t k>) (4 (t k> - £k<V) 

+ 0(fi)]/24 (All-25) 

Introducing (A7-18) and (A7-19) gives 

-ko = -k (t " V “ ^k (t ‘ V / 2 + [ G k -k 

-*c] (t - v 3 / 6 + - PkM] - v 4 / 24 

+ o((t - t k ) 5 ) + o(p(t - t k ) 4 ] (All-26) 

(All-26) represents the behavior, in the overlap domain, of the zeroth order 

til 

relative position between the trajectory and the k body. It is used to 
determine the behavior of the higher order terms and in the matching. 

All. 2 First Order 
Let 

r, = X + X. + X_ + X- 

— ko — o —1 —2 —3 

where 


So = V - ‘k> 

(All-28) 

Xj = -K£ k (t - t k )/2 

(All-29) 

*2 = [G k V k - <t - t k ) 3 /6 

(All-30) 

—3 = 0 ((‘ - l k> 4 ) 

(All-31) 
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Then, using (A2-14) 


t(£ko> = t<2o> + G <*o> [*l + —2 + °<*3>] 

+ HfX^ [xj + 0(X 1 X 2 )]/2 + O(Xj) 

In order to use the solution being developed for both approach to and 
from the k** 1 body, let 

(t - t k ) 

Q k = sgn (t - t k ) = -p— j-j 
T hen 

i<* o> = i(V 1 - < k >) 

= -I k (‘-t k )/|v k (t-t k )| 3 

= Q k KYfct/lt - t k ) 2 

By a similar analysis using (A2-15) and (A2-16) 

°<x,) = Q k c(Yk>/“ - V 3 

S<*o> = Q k H(V k )/(t - t k ) 4 

Now (All-33) can be written 

l(£ ko ) = —7-2 [K^t - I G <V k > £ k (*'- * k > 

(t - t k ; 

+ |G(V k ) (G k V k -„ £k : ) (t-. k ) 2 
+ °(<‘ - ‘k' 3 )] 


(All-33) 
de parture 

(All-34) 


(All-35) 

(All-36) 

(All-37) 
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(All-38) 



The first order force function is given by (A5-9), It can be rewritten in 
terms of a singular part, F^ g , and a non-singular part, Fj n , i.e., 

El(£ 0 . £i> = F ls (V 2k 1 + tin'V £i> (A 11-39) 

where 

£ k ) = M k l(r ko ) (A 1 1-40) 

N-2 

F ln (£o. £i> = M k tJSk> + 2 - £l> + i<£i>] (All-41) 

1 = 1 
i^k 

Using (All-38) and (All-40) the singular force function is 
^k f u * 

lis = 777[l( v - k )-fG(V- k )£ k ( t-t k ) 

u - t k ) 

+ 1 G(V k ) (G k V k - p£**) (t - t k ) 2 + o((t - t k ) 3 )] (All-42) 

From (A3-28) 

B (t, t) = I(t - .) +G(t) (t - T) 3 /^ + o((t - t) 4 ) (All-43) 

whe re 

G (t ) = G^t)) (A 11-44) 

According to (A6-11) the first order solution is a function of B(t, t) Fj(t). 
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From (A 11 -42) and (A 11 -43) 


B<t, T) F ls (T) = fi lfc -^^2 + £ 2k + —3k <* - T » 

(T " V ^ 


+ p - + o 

+ —4 k . t ,2 + u 

( T - l k ) 


(t - t) (t - t k ), 


4 3 

(t - t) (t - rf 

(T -t k ) 2 ’ "« T -V 


(All-45) 


whe re 

£ lk = Q k X(V k ) (All-46) 

P 2k = ^ Q k M k G <^k' E k / 2 < AU - 47 » 

P 3k = Q k M k G( V (G k - ^V 6 (All-48) 

—4k = $4k (t) 

= Q k M k G(t) £_(Y k )/6 (All-49) 

(All-45) represents the singular behavior of the integrand in (A6-11). 

Although only the first term of (A 11 -45) is actually singular as t -*t^ and 
t -*• t k the other terms play an important role in evaluating the second order 
solution. It is not necessary to develop a corresponding expansion of the 
non- singular part of the integrand in (A6- 1 1 ), i. e. , B(t, t ) .F ^(t), as this 
contribution can be evaluated directly as will be shown. 
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Now let 


*Ik (t ’ T) = B(t - T) F -ls^ - filk - fc=k (TTTJ 

(t- t k ) k' 


P 3 k (t - T) “ hk {t) 


(t - t)~ 

(t_ v 2 


(All-50) 


From (All-45) 


± lk (t, t) =. O 


(t - t) (r- t ), (t ‘ J ^~2> M. £ “ t T) } 

( t -V { k 


(All-51) 


The singular part of (A 6-11) is 


/ 


B (t, r) F ls (t) dr = 


t 

{ 


£lk ( " ■ ja + 2.2k (t- t. ) + 2.3k (t " T) 

(t - t fc ) v k 


+ + iik ^ 

( t - t. ) 


dT 


(All-52) 


At this point (All-52) is exact since the last term on the right integrates the 
error in the expansion (A 1 1-45). The first four terms on the right integrate 
to 



, (t ' t, .2 dT = - >°« Q k <‘ - *k> + >°s Q k (‘ 0 

(t- t k ) 



- 1+(t - t k , / (, o- t k> 


(All-53) 
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[ (r-"t J dT " (t o ‘ t k ) + (t ' t k ) l0g Q k (t V 

t k 

o 

-j^+loB Q k (t 0 - t k ) ] (« D - t k ) 

/ <* - -I dx = (t o - .//> - (t o - t k ) (t - t k > 

'o ' 

+ (t - *k>7 2 

A 3 

/ d ' = < l o - V 2 / 2 - 3(t o - «k> <* - ‘k> 

t (T " t k ) 
o 

- 3(t - t k ) 2 log Q k (t - t k ) 


(All-54) 


(All-55) 


+ 3 [l/ 2 + log Q k (t 0 -t k )] (t-,/ 

+ (t - t k ) 3 /(t o - t k ) (All-56) 

Note that (A 11- 53) contains a logarithmic singularity as t -* t k> This is the 
singular behavior of (A 6-11). On the other hand, the first two terms in 

g 

(All-50) cancel out the singular behavior of B(t,T) F. (t) and is finite 

— iS 1 K. 

as t t k> Therefore the last term in (All-52) can be expanded about t = t k> 
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The form of the expansion is 


f 


-lk T ^ dT 


f ilk (t ' T 


)d t 


t=t. 


5 / ilk <*• 


t)c3t 


( t - t k ) 


t=t. 


- t 

at 2 J 


$ lk T ) dr 


(t - t k ) / 2 


t=t, 


+ O 


(“ - l k> 3 ) 


(All-57) 


The individual integrals in (A 11- 57) can be evaluated as follows: 




(t, T)di 


J t=t, 


/ 


ilk (t k' T,dT 


/ k b 


T) F. (t) + p. 


-Is^' "^lk (T- t k ) 


dT 


&./ 


+ Pot. / dT + ^3k+ / ( T - t k^ dT 
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+ P 4k (t k } / (T ' t k )dT 
*0 

r k s 

= J — 10k (t k* T)dT " —2k (t o " V 


- [>3k + i4k <V] <‘o - '/A (A 11 -58) 

where 


ilOk <V T > = B(t k' T >t la M +£lk/(T- t k ) (All-59) 



■V 

J 


1 

V 

_l 

d 

dt 

/ i“k <*• T,dT 

II 

J silk (t ’ T)dT + iik (t > 


•/ 

t 


i 

- 

O 

t=t, 

o _ 


/ 


^ii k (t ' T ) dT 


(All-60) 


since, by (A 1 1 -5 1 ) 


£ lk (t, t) = O 


(All-61) 


Using (A3-16) and (All-50) 


di 


lk 


dt 


= D(t, t) F 1s (t) - 


£ik 


(T - t k ) 


_±2k_ _ 

(r- t k ) —3k 


- 3jP 4k (t ) ~ ——————— — 


( t - t,)‘ 


dt {t ~\) Z 


(All-62) 
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f tk <3 ,■ „ H 

J dt ^lk (t k’ T)dT = 

t 

o 



D(t k . t)F u (t) - 


— lk 

< T - V 2 


■t-2k 


{t - t k ) J 


dT 



o 


+ 3 £4k<‘k» - i &4kW T - 'k> 


dT 


/ ** 

*!ik (V T)dT + [p 3 k + 3 ^k <‘k>] (* 0 - V 

t 

o 


- Jf^kMV l // 2 


(All-63) 


whe re 

*Uk <V T) = D( V T » £u (T) - £lk /,T - 'k> 2 - £ Zk /(T - V (All-64) 

Finally 


2 r* 

/ ll k <t.T)dT 

dt -/ 

II 

f TTi?k (t * T)dT + J^ik (t ’ t} 

•J dt 

( 

1 

r+ 

0 

1 

t=t,_ 

r+ 

0 

1 



t 

o 


7J iik (t k ,T)dT + ° <|X) 


(All-65) 
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since, by (All-51) 


■St T > 


o 


(T - V 


(t - T) 

( T - t k ) 


3 

2 * 


(t -t) 
(T - t k 


^*lk (t > -r) = ©[(t - t k ) , (t - t k ), jx (t - t k )] 


it iik (t k’ t k ) = ° (t ■ t k ) 


(All-66) 


Using (A3 - 18) and (All-62) 


j 2 s 

d ilk 

dt 2 


G(t) B(t, r) F lg (r) - 6P 4k 


(t - T) 
(T - t k )‘ 


dP 


— 4 k (t - t ) 2 d -^ 4 k (t - t ) 3 

dt (T - t k ) 2 dt 2 (T - t k ) 2 


(All-67) 


J dt 2 


lk (t k' T) dT = 


n 


P4k (t k ) 


G (t k ) B (t k , T ) F 1 s (t) + 6 (T _ t ) 


dT 


■/ 


6 dt 14k <*k> - f? -4k<*k> < T - ‘k> 

dt 


dT 


f l k 

J iuk<V T,d + 6 JfJW‘ 


k» <*o - V 


-fj£4k<V (, o- V 2 / 2 

dt 


(All-68) 
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where 


il2k (t k’ T) = G ^ t k ) B T) B 1s { t) + 6 -^k <V/< T - V (All-69) 

Using (A 11-49) and then (A 11 -44) and (All-2 5) 

-12k (t k’ t) = G k— 10k (t k’ t) (All-70) 


Substituting the results of (A 1 1- 58) through (All-70) into (All-57) gives 


-lk dT = 


-10k (t k- T > dT ‘ £2k <*o ' *k> " fck 




+ [—3k + 3|3 4k (t k’] <*„ " V 

r tk ■ , 

+ G k J -10k’ (t k’ t) dT + 6 dt -^4k (t k )(t o " V 

t 

o 

'-fr^kMV-'k ' 2 / 2 <‘-V 2 / 2 

dt 

+ o((t - t k ) 3 ) (All-71) 
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The derivatives of i^(t) are simply coefficients in the expansion 


£ 4k (t) - + dt -^4k ( V (t " V + j 2 £4k (t k^ (t " V / 2 

at 


+ o((t - t k ) 3 ) 


(All-72) 


Substituting (All-53), (All-54), (All-55), (All-56), (All-71) and (All-72) 
into (A 11- 52) and collecting terms gives 

/ T > « dT - - £ lk 0 k « - V *£| k [ lo s Q k <' n - V - 1 ) 


r 


+ 1 — 10k (t k* T) dT + —2k (t - t k ) log Q k (t " t k ) 


£lk/ (t o " V " £ 2 k[ log Q k (t o ' V + *] 


/ 


ink (t k' T)dT 


(t - V 


3 - 64 k ' V log Q k <( ' V 


—3k + 3 &4k ^b 1 + 6 ^4k (t k’ log Q k (t o ' V 


ft 


+ G k I iiok ( V T)dT 


(t . t k r/2 


+ o 


( (t - ‘k' 3 ) 


(All-73) 
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(All-73) is an expansion of the singular part of (A6-11). Expansions of the 
non-singular terms in (A6-11) can be obtained using some of the expressions 
found in Section A3 and Taylor series formulas.. The results are 


A(t, t Q ) = A(t k . t o ) + C(t k . t Q )(t-t k ) 


+ G k A(t k , t Q ) (t - t R ) jz. + o((t - t fe ) ) 


(All-74) 


B(t, t Q ) = B(t k , t Q ) + D(t k> t Q ) (t - t k ) 


+ G k B(t k , t o ) (t - t k ) 2 /2 + o((t - t k ) 3 ) 


(All-75) 


f B(t, T) £]„(T)d T = f B(t k , T) F ln (T)dT 


/ 


D(t k , T) F ln (T)dT 


(* - t k ) 


/ 


G k I B(, k , T )F ln (T,dT 


whe re 


+ — ltAk^ 


(t - tj.f/2 + o((t - t k ) 3 ) (All-76) 


— ln (t k } -?k 


(All-77) 
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Substituting (A 1 1-73) through (All-76) into (A 6.- 1 D.gives 


£_1 (t) = 

a lk log Q k (t - t k ) + b lk + lk (t - t k ) log Q k (t - t k ) 



+ ii k (t ' V + lk (t " t k ) + ilk (t " V log °k (t " V 



+ i-lk (t " l k )2 + °( (t ' V ) 

(All-78) 

whe re 



-lk 

p „■ L 

= “ilk 

(All-79) 

^lk 

= i lk [ lQ g Q k (t o ' V _1 ] + A ( V t o ) ^1 <V 



+ B (t fc , t Q ) vj (t Q ) + K 10k (t k , t Q ) 

(All-80) 

^-lk 

= i2k 

(A 11-81) 

iik 

= Pi k /(t 0 ri V + C V V-l (t o ) + D (t k’ t o ) il {t o } 



+ Kllk<v t G ) 

(All-82) 

lk 


(All-83) 

ilk 

• - 3 P4k (t k» 

(All-84) 

&lk 

= P 3k /2 + 3 £ 4k <‘k> [ log Q k <‘o - V + 1/Z ] + £ k/ 2 



+ G k [ A (t k . t 0 ) Li (to) + B (t k’ t o ) il (t o } + -10k (t k’ 



(A 11 -85) 


{ \ 
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and 


-10k (t k’ V " j [-10k (t k’ t) + B (t k’ t) -In (r) ] dT (A 11 -86) 

t 

o 

Silk (V V = f k [*!lk <'k' T > + D <‘k’ T > —In ' T )] dT (AU - 87> 

*t 

o 

(A 11- 78) gives the behavior of £^(t) in the overlap domain. The logarithmic 
singularity as t -*• t k produces the non-uniformity in the first order outer 
solution mentioned at the end of Section A6. It will be shown that a similar 
singularity occurs in the inner solution. The integrals given by (All-86) and 
(All-87) cannot be evaluated analytically, i. e. , expressed in closed form. 
They must be evaluated numerically and this problem is discussed in 
Section C. 

All. 3 Second Order 

In deriving the first order expansion it was first necessary to derive f (_r ko ), 
c. f. (All-38). According to (A5-10) the second order force function requires 
an expansion of G (r ^ ) . From (A2-20) and (All-27) 

G (£ ko ) = G ^o) +H(x o )[x 1 +x 2 + 0(x 3 )J 

+ T (£o>[*i + 0(xj x 2 )]/2 + O(x^) (All-88) 

Continuing the sequence (All-35), (All-36) and (All-37) gives 

T^) = Q k T (V k )/(t - t k ) 5 (All-89) 

‘I 
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Then (All-88) becomes 


G t k o> 


*> - 2 2 <V k ) £ k (t - t k ) 


+ l 5 (V k ) (O k V k - «-**) (t - t k ) 2 


+ £ I<V k )£ k 2 («-t k ) 2 +o((t-t k ) : 


(All-90) 


The second order force function can be written in terms of a singular part, 
P\, , and a non- singular part, F^ n , i. e. , 


I 2 (£o. £ r £i> = £ 2s (£<,■ I,. £ k > + Z Zn (Jo- Pi) 


(All-91) 


where 


II- Ek> = M k G <£ ko )£i 


(All-92) 


F (r , r., p.) 
— 2n — o —1 l 




IN -C 

L\ + 2 M i G { - “£i>£i (All-93) 


Since is a non-zero vector, Hjr^) is non- singular. The dyadic _r^ 
contributes a factor of (log (t - t^)) 2 but 


[(log x) 2 dx = x (log x) 2 - 2x log x + 2x -»• O 


as x — O. Therefore the vector H(r ) r. is finite as t -*-t, and is included in 

, — — o —1 k 


F^ n (also see Carlson, pg. 27). 
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Multiplying (All-78) by (All-90) and substituting into (A 11- 92) gives 


* lo § Q k (t ~ t k ) G k -lk 

—2s _ G k — lk 3 + 3 

2 k 1 (t - t k ) 3 (t - t k ) 3 


+ M- (G‘ c„. - J,‘ a„J 


log Q k (t - t fe ) 


k^lk ' “k-lk' . .2 

1 ' V 


+ 


G k -lk + ^ (G k -lk - J k 


- lk> ] <t -\j- 


G k -lk + H k -lk ' ^ J 


2 * 1 lo e Q k (t ' t k ) 


k -lk 


(t - t k ) 


2 * 


G k ilk + H k -lk " * J k -lk " ^ J k -1 


lk 


(t - t k ) 


+ 0(log It - t. 


where 


G k = M k Q k G ^k> 


H, 


M k Q k 


IW ( G kVk-^VS W£k 2 


M k Q k 


J k. = -V^^Pk 


According to (A6-12) the second order solution is a function of (B(t, i 


(All-94) 

(All-95) 

(All-96) 

(All-97) 
) F 2 (t). 


52 



From (All-43) and (All-94) 


Bit, r) F 2s (t) = i lk i°g Q k (f- t fc ) +i 2k 


< T - ‘k r 


( T - 'k>' 


where 


(t - t) 


(t - t) 


+ ^3kT i ~ L J l0g Q k (r " t k ) + —4k , . ,2 


< T - V 


<T - ‘k 1 


+ 4 5 kf J T7 log °k (T - ‘ k > + 46k ^77 

( T -V (T -V 


3 3 

+ +7k iL ^ L 3 log Q k<‘ - V + 4gk 


(T - V" 


(t - t k r 


+ o 


(t - t) log |T- t k | , ■ |t ~ - 1 log | T - t k |, 

• lr - V 


(t - t) 
(T - *k r 


l0g | T - *k 


(All-98) 


^lk = G k-lk 


(All-99) 


—2k = G k^lk 


—3k 


( G k £i k ' J k £ lk ) 


^4k ' G k -lk + ^ (G k -lk ' J k -lk ) 


—5k G k -lk + H k -lk " ^ -it' 


k -Ik' 


(All-100) 
(All- 101) 
(All-102) 
(All-103) 
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(All-104) 


- 6 k = G k —lk + H k -lk ' ^ J k -lk " 11 J k 2.1k 

-7k = —7k (t) 

= 0(t) G*a lk /6 

isk = t 8k (t| 


(All-105) 


(All-106) 


(All-98) represents the singular behavior of the integrand in (A6-12). It is 
not necessary to develop a corresponding expansion of the non- singular part 
of the integrand in (A6-12), i. e. , B(t, t) F^ n ( T )> as this contribution can be 
evaluated directly. 


Now let 


— 2 k T » 


(t - t) 

B (t, T) F 2s (t) - £ -5 log Q k (T - t k ) 

( T - t.) 


—2k - —3k log Q k < T - V 

(T -V {T -»k> 


± 4 k - ** ■ r , T \i - i 5 k 777 106 Q k (T - ‘k> 

( T - t k ) (T-t k ) 


— 6 k - ±7k <‘> — 7^ log Qk < T - *k> 

k (t - t k ) 


i S k<‘>f7% 

< T - l k> 


(All-107) 
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From (All-98) 


tzk <*• T) = ° 


(t-T) lOg | T - t fc 


(t - T) 


(T- t k ) 


j!og | T - t k 


(t - r) 

< T - 'k r 


log 


(All-108) 


The exact expression for the singular part of (A6-12) can be written 



B (t> t) F 2 s (T)dr 



JLiU' iogQ k ( T -t k ) + i 2k 

(T - fc k } 


(t~T) 

( T - t k ) 3 


+ i,,. ■ (t - ~ ~2 log Qt (T - t J + * -<*- r) - 


—3k , , ,2 ~k ' ' ~ -k 7 •' —4k , , ,7 

( T - t.) (r- t, ) 


+ 4^ 7_ - Vr iog Q, (t- t. ) + 4?,- 


-5k (t- t k ) 


k ' ' -k 7 1 —6k (x - t fc ) 


+ “7k ,t) 7 1 ~^ logQ k (T - t k ) 

(T - 'k> 


+ —8k (t) 


(t - T) J 

(T - ‘k> ; J 


c3t 


+ 


/ ilk <*• T) dT 

t 

o 


(All-109) 
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The first eight terms on the right of (A 11- 109) integrate to 


/ 


(t - t) 




3 log Q, (T - t k ) dT 


[log Q k (t - t k )] 


2(t - t k ) T 4(t - t k ) 


N Q k <‘o - V + *] 


(t - t, ) 
o k 


[ 2 lo 8 Q k <*o - V + *] 


4(t - t, ) 
' o k' 


2 (t ' V 


(All-110) 


f . I* - T ) dT ! I + . 

I <-**>’ 2(t 'V (t o-V 2 <t 0 -V 2 

o 


(All-111) 


/ 7^77?! log Q k (T - y dT = 


. <T -*k> 

o 


- J log 2 Q k (t - t k ) - log Q k (t - t k ) 


+ j log 2 Q k (t Q - t k ) 


[log Q k (‘ - ‘ k > + l] 

1 + t - t V 

1 o k 


(All- 112) 


/ 


(t - t) 


i (T - *k> 


2 dT = - lo « Q k (t - *k> + log Q k (t o ‘ V 


- 1 + 


ft- V 

(t o - V 


(All-113) 
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/ rr^T) 106 Q k (T - V dT = (t o 1 V [ log Q k (t o - V - '] 

t k 

O - 

+ j (t - t k ) log 2 Q k (t - t k ) 

- {t - t k ) log Q k (t - t k ) 

- [jlog 2 Q k ('o ' l k > - 1 ] ('-V 

(All-114) 




t 

o 


dT = «*o ' V + <* - V log Q k <• - V 

- (log Q k (t Q - t k ) + l](t - t k ) (All-1151 

3 * 

, (t ' \ log Q k <T - l k> dT = (t o - V [ log Q k (t o - V - M 

( T - t k ) 

+ § (t - t k ) log 2 Q k (t - t k ) 

+ § (t - t k ) log Q k (t - t k ) 

- [I *° gZ Q k (t o - *k> - t] < l - V 

+ o((t-t k ) 2 ) (Ali- 116) 
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(A 11- 117) 



(t - T)° 
(T - 


dT 


(t o - V + 3 ft ■ V l0g Q k (t ' V 


-[3 logQ k (t 0 -t k )-l](t-, k ) + o((t-t k ) 2 ) 


Note that (A 11- 110), (A 11- 1.11), (All- 112) and (A 11- 113) contain singularities 
as t t^. This is the singular behavior of (A6-12). On the other hand, the 
first five terms on the right side of (A 11- 107) cancel out the singular behavior 
of B (t, t) F ^ (t) and $® k is finite as t 1^.. Therefore the last term in 
(All- 109) can be expanded about t = t^. The form of the expansion is 


r* 

Y 1 

J £2k (t * T)dT = 

/ -2 k ft’ T)dT 

t 

0 

J 

t 

^ 0 J 




$2k ft* T ) dr 


t=t, 


ft - t k ) 


+ O 





(All-118) 


The individual integrals in (All-118) can be evaluated as follows: 



A 

J —2k <*k' T)dT 



t 

o 


B(t k , t) Y _ 2 s (t) + vp lk 


log Q k (t - t k ) 

< t - t k > 2 
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1 


+ 4 > 




log Q k (r - t k ) 


^ {t _ v 2 ~3k (t - t k ) 


+ 4- 


“4k (t - t k ) 


dT 


+ [+5 k + —7k <' k >] / l0 S <\ < T - ‘ k > * 

t 

o 

+ [46k + 4g k <‘k>] / dT 

t 

o 

= / *20k <V T) dT - [*5k + i?k <‘k>] "o - ‘k> 


[l°g Q k (t c - t k ) -l] - [+ 6k + 4g k (t k )](t D - t k > 

(All-119) 


whe re 


log Q k (t - t k ) 

iaok ('k- T » = B,t k' T > -2 s ( t) + 4lk ~ t ,2 - 

<T - V 


l lo § Q, (T- t k ) , 

+ 4^2 k 2 + ^3k Tr — tl + ^4k 

(T-t k ) Z 3k (T -V 4 k ( T -t k ) 


(All-120) 
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r t 


A 

dt J —2k (t ’ T) dT 

t 

0 


. S -2k (t ' T) dT + ilk (t - *> 

t 

O 


= / ^^2k (t ' T)dT 

t 

o 


(All-121) 


since, by (All- 108) 


i 2k (t. t) = O 


(All-122) 


Using. (A3-16) and (All-107) 


d s logQ (r_t ) 

dt ^ 2 k <«’ T) = D(t. r) F 2 s (t) ■ », k — ~ ~ 3 ~ 

(t - t k ) 


- 4, 


1 , l0g Q k (T - 




—2k 3 — 3 k 2 

( T -t k ) (T-t k ) 


, 1 , lo 8 O k ( T - t k ) 

^k-^-^z-isk (T-t^ 


- +6k 7—7- - 3 *7k <‘> 7-77 «k (T - V 

( T -t k ) - (T-t k ) 


3 -8k (t */ I 3 _ dt —7k ^ ^ ^ 3 log Q k (T - l k' 

(t- t k ) (T- t k ) 


Ad, (t ) (* - m 

dt -t8k l , .3 

(t- t k > 


(All- 123) 
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f tk si 

J * 


A. s s 

dt -2k (t k’ T ) dr = 


/ 


* 2 i k <v T > dT -|-i 7k (t k )(, o . tk) 




where 


(All-124) 


*2 Ik ( V T > = dt -2 k (t k ' T > -1^) log Q k (T . t } 


" dt -8k { V 
Substituting (All_ H9) 


(All-125) 


and (A 11 - 124) into (All_n 8)gi 


gives 


A 

J — 2 k T ) dT = f $ s , r 

t J - 20 k (t k' T > dT - U 

O t L- 


5k + —7k (t k } J (t o - V* 


[ logQ kV- 


/ 


5za ( ^ T,dT 4 i 7 k<V(*o-v 


!°g Q, (t - t ) ,1 d 
L k 0 k rdtiskty^-y 

+ («-v 2 ) 


<t - t k ) 


(All-126) 
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The derivatives of and are coefficients in the expansions 


^7k (t) -■ —7k ( V + dt ^7k (t k )(t “ t k^ + ° l (t " V ) 


(All- 127) 


-8k (t) —8k ( V + dF-8k (t k )(t “ V + °( (t V ) 


(All-128) 


Substituting (All-110) through (All-117) and (All-126) through (All-128) 
into (A 11- 109) and collecting terms gives 


f‘„,. , ±lk lo 8 Q k (t -‘k> . 1 Phk .4,,,) 1 

J B(t> T) —2s (t| dT ' 2 (t - t ) + 2 \ 2 + 2k / (t - t. 


-^log 2 Q k (t - t k ) - (+ 3k + + 4k ) log Q k (t - t k ) 


- + 


[log Q k (t„ - t k ) + l] +. 


2 k 


lk (t - t, ) (t - t. ) 

ok 'ok' 


+ +3k 


I log 2 Q k (t D - t k ) 


- 1 


+ —4k 


log Q k (t o - *k> - 1 




20k (, k’ T > dT 


1 

+ 2 


fek + 3 ±7k <‘k> 


(t - t, ) log Q, (t - t ) 


—5k - *6k - 2 —7k <V - 3 *8k <V 


(t - t, ) log Q, (t - t ) 
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iik 


[2 log Q k (t o - t k ) + 1 ] +. 


2 k 


4(t " t, ) 
o k 


2(t o - ‘k> 


+ i3k 


[iog Q k (t c - t k ) + l] 4^ 

<‘o - *k> ETV 


- +5k.[j log2 Q k <‘o - *k> - *] 

- itk [ 1o « Q k I'o - 'k 1 + '] 


[I lo s‘ 


^7k (t k) Lf l0 8 2 Q k <*o - V ' X 


‘k ' o V 4. 


" -8k ( V [ 3 log Q k (t o " t k ) " *] 


fa 

J —2 lk {t k’ t) dT 


( t - t k ) 


+ O 


( {t - V 2 ) 


(All-129 


Apart frorrybeing the longest expression derived thus far, (A 11- 129) 
represents the singular part of (A6-12). Expansions of the non- singular 
terms in (A6-12) can be obtained using some of the expressions found in 
Section A3 and Taylor Series formulas. The results are given by (All-74) 
and (A11 -t 75) as well as 


/ 


B(t, t) F, (t) dr = 


/ 


B (t k’ T) — 2n (t) dT 


/ 


D(t k’ T) — 2n (r) dT 


(t - t k ) + o((t - t R ) 2 j 

(All-130) 
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which is just (All-76) with replaced by F . 

Substituting (All-74), (All-75), (All-129) and (A 1 1 - 130) into (A6-12) 
gives 

—2 ^ = —2k ^ - ^k ^ ”• *k^ ^ —2k ^ ~ ^k^ 

+ £ 2 k log Q k (t - ‘k> + —2k 

+ e 2k (t - t k ) log 2 Q k (t - t k ) 

+ L 2k <* - *k> log Q k <* ' V + —2k (t - V 

+ o[(t - t k ) 2 , ^log 2 |t - t k |] (All-131) 

where 


a 2k = + lk /2 (All-132) 

i 2 k = < 3 llk + 2 i 2 k>/ 4 (All-133) 

c 2k = -+ 4k .(All-134) 

—2k = - ilk [ log Q k (t o - V + *] < l o - V"' + i 2k (t o - V" ‘ 


+ Aik [ log Q k <*o " t k ) ‘ ! ] + -20k (t k' ‘o* 

+ A ^k’ t o^ -2 (t o 5 + B ^k’ t o^ -2 (All-135) 

A>k * [+5k + 3 *7k < l k>]/ 2 (AU - l36 > 
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(All-137) 


—2k = “ [ 2 -5k ' 2 -6k " 3 — 7k ' 6 — 8k (t k } ] / 2 
i2k = ^lk [ 2 log Q k (t o " V + X ] {t o " V" 2 / 4 
+ ^2k (t o " t k ) /2 + i4k (t o ' V 

- i 5 k [ log2 " t k ) - 2 1 / 2 

- i6k [ log Q k (t o - V + 1 ] 

~ 7k <V [ 2 l0g2 Q k (t o - V - 5 ]/ 4 

- i8k (t k } [ 3 log Q k (t o " V " *] + -21k (t k’ t o ) 

+ C ^k’ t o^ -2 (t o ) + ° ^k’ ^ -2 ( V (All- 138) 

and 

K 20 k (t r’ ‘o' = / [*20k <V T > + B (t k' T > < T >_ * (AH-139) 

—21k ^k’ t o ) = j [-2 lk (t k’ + D (t k’ t) — 2n dT (All- 140) 

(All- 131) gives the behavior of r^t) the overlap domain. Like the first 
order solution it contains singularities as t-*t^ which must be cancelled by 
similar singularities in the inner solution. Also, just as in the first order 
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expansion, the integrals (All- 139) and (A 11- 140) must be evaluated 
numerically; a problem which is discussed in Section C. 

All. 4 Third Order 

The third order solution is not actually used but knowing its behavior in the 
overlap domain is important in the matching. From (A5-8) 


£3 = G( £o^ 3 + —3 — 1’ —2’ £i> 


(All-141) 


with a solution similar to (A6-11) and (A 6 -12) 


£3 (t) = A (t, t Q ) £ 3 (t Q ) + B (t, t Q ) v 3 (t Q ) 


/' 


+ I B(t, t) F 3 (t) dr 


(All-142) 


Carlson has shown that in the overlap domain 


-3 


■(’ 


r, (t) = 0((t - t k )' 2 log 2 |t - t k |j 


(All-143) 


A 12 BEHAVIOR OF THE INNER SOLUTION IN THE OVERLAP DOMAIN 


A12. 1 Zeroth Order 

According to (A9-1), (A9-3) and (A9-4) the zeroth order inner solution is a 
function of the eccentric anomaly F^. The behavior of the outer solution 
in the overlap domain has been developed as a function of time. Therefore 
it is necessary to find the behavior of the inner solution as a function of 
time, i. e. , find 


F 


k 


F k (S k» 


(A 12- 1 j 
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Such a result is found by inverting (A9-2) but, in general, this is not 
possible due to the transcendental nature of the equation. However, the 
overlap domain for the inner solution corresponds to the region far out on 
the asymptotes of the zeroth order hyperbola. This region is characterized 
by large values of the time and this makes it possible to invert (A9-2). 

In order to be completely general it will be assumed that on the approach 
asymptote of the inner hyperbola and F^ are negative and on the departure 
asymptote they are positive. This will require two solutions of (A9-2). 


Dropping the subscript k for the time being, (A9-2) can be written 


nS = (ee F - ee' F - 2F)/2 (A12-2) 

\ 

For F— *oo (A12-2) is rewritten as 


e F = (2n S + 2F + ee F )/e (A12-3) 


_ J* £ 

Since F and e are much smaller than e (A 12-3) can be approximated by 


e F = 2n S/e + O(F) (A12-4) 


and then 


F = 



(A12-5) 



(A 12-6) 
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Substituting (A12-5) and (A12-6) into (A12-3) gives 


e F * l|S + | log ( 2 g§)+ o fiats] 


(A 1 2-7) 


(A12-7) is simply a more accurate form of (A12-4). From (A12-7) 


F - log (#) + i log (^) + „ 


(A 1 2-8) 


r F = + o ( L°-g,s 

** U 2 / 


(A 1 2- 9) 


Substituting (A 12-8) and (A 12-9) back into (A 12-3) gives 


2nS , 2 , ZnS\ . 2 , / 2nS \ „ log S 

— +=log (— ) +==s l °g (— )+**? + 0 (~V“ 


2nS 


(A12-10) 


F = l og ( 2|S)+^l og 



2 2 
4n S 


2 log' 



. , / 2nS \ -2 

41og ( — - e 


+ 0 


' 2 

log 5 

c 3. 


(A 1 2 - 1 3 ) 


-F 


e 

2nS 


2n S 


. / 2nS 

2„2 log “ 


+ 0 


' 2 ' 


(A 12-14) 


(A12-13)is the explicit form of (A12-1). Two additional useful expressions are 


sinh F = (n S + F)/e 


cosh F = sinh F + e 


-F 


(A 1 2- 15 ) 
(A12- 16) 
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For F — -oo, (A12-2) is rewritten as 


e" F = - (2nS + 2F - ee F ) / e 


(A 1 2- 17 ) 


The same sequence of approximations as used for F positive leads to 


-F 


2nS . 2 


2ns\ 2 


2nS 




log I si 
, s 2 . 


(A 1 2- 1 8 ) 


, 2nS \ . 1 , / 2nS \ . 1 

F - lo e I-— ) + ls l °S (- —j + -^2 


2 log' 


2nS 


- 4 log 


r 2nS \ -2 


+ 0 


log I S | 

^3 


(A12-19) 


-F 


2rTS 


2~2 log 
2n S 


2nS 


+ 0 


- lQ gl S l 

r.3 


(A 12-20) 


The function sinhF is still defined by (A 12- 15) but coshF can now be written 


coshF = - sinhF + e 


(A 1 2- 21 ) 


The departure asymptote with F > 0 must overlap with the outer solution where 
t - t^ > 0. Likewise, the approach asymptote with F < 0 must overlap with the 
outer solution where t - t^ < 0. Introducing defined by (All-34) makes it 
possible to write one expression for F: 


2Q,n,S, \ / 2Q, n S 

= log ( 1"^) + — log -^ k 

k / n k k \ k 


Q, 


—2 2 


21og2 (!2LAV-ogf^)-^ 


+ o 


log 2 |S, 


(A 1 2- 22) 
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Using (A12-22) in (A12-15) gives 


sinhF, 


^ + £k lo 

e k e k ° s \ e k 


A 1 . I “"Vk k 

n, e, S, ° g 1 e, 


Q, 


4-2“ c 2 

4n k e k S k 


2 log 


2/ 2Q k n k S k 


/ 2 ^k n k S k \ —2 

4log( — V— 1- e ? 


+ 0| 


log 2 | S, 


(A 12- 23) 


Using (A12-14) in (A12-15) and (A12-20) in (A12-21) gives 


c ° 3hF k ' Q k sinhF k + 


k“k Znffi 


. 2Q n S 

k , / k k k 

log 


+ 0 


log 2 | S, 


(A 12- 24) 


(A12-22) through (A12-24) represents the first step in the development of the 
behavior of the inner solution in the overlap domain. 

I I 

Now consider a rectangular x, y coordinate system in the inner orbital plane 

I 

defined by If the x axis lies along the line of nodes between the orbital 

plane and the x, y reference plane with x positive toward the ascending node 

i ... . 8 

then the motion is given by 


— 1 i A 1 R'ri 

x = A - — coshF - — = — sinhF 
e e 


(A 12- 25) 
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(A 12-26) 


ni a ' C 1 

y 1 = B' - — coshF + — : — sinhF 
1 e e 


where 


A 1 = e a cos go 


B' 


e a sm to 


C' = (e 2 - 1) 1/2 


(A 12- 27) 
(A 12-28) 
(A 12-29) 


Here to is the argument of pericenter measured in the orbital plane from the 
ascending node. The subscript k has again been temporarily dropped. 

Using (A 12-24) x' and y' become 


x 


Q 


' - - (A* + QB'C') sinhF 

e 


+ A 1 


1 " 27TS + ,-2 2 log 1 
2n S 


1 2Qns\ / log 2 ISp 

^ / \ i 


(A 12-30) 


y' = - % (B 1 - Q A ' G 1 ) sinhF 


+ B' 


Q ■ ^ 


i - — =— + 

2nS 2„2 
2n S 


7 


(A 12-31) 


Differentiating with respect to S gives 


— = -■= (A 1 + QB'C') coshF ~ + 0[-U 
dS e dS \s 2 / 


(A 1 2- 32) 


4X- = _ -8- fR' - 


, Q - --^(B 1 - QA'C') coshF 0(4y 

dS e dS Vs / 


(A 12-33) 
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Differentiating (A9-2) gives 


— dS _ — , „ , 

n -7T7 = e coshF - 1 
dF 


(A 12-34) 


or 


dF 

dS 


n 


e coshF - 1 


(A12-35) 


Then 


coshF 


dF n coshF 


dS e coshF - 1 


(A 12-36) 


As F — ► co, S — ► oo and 


i -rp dF n Q sinhF 

coshF -jrr — — 


dS e Q sinhF - 1 


n 

"e 


(A12-37) 


because of (AT 2- 23). Therefore the velocity components as S — ►co'are given by 


u' 

u dS 


= - Q n (A 1 + QB'C 1 ) / e‘ 


n (B'C’ + QA') / e 2 


(A 12- 38) 


V = 


dy ' 
dS “ 


Q n (B' - QA'C' ) / e 


-2 


n (A'C - QB') / e 2 


(A 1 2-39) 
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The components U 1 and V' are the orbital plane components of the hyperbolic 
excess velocity. V . 

CO 

(A 12- 30) and (A12-31) can now be written 


— U'e 

x' = sinhF + A' 

n 


Q 

2nS 


+ -2 2 log 




(A 12-40) 


— V'e 

y' = ■ sinhF + B' 

1 n 


-_Q_ + _ L_ 

2n S 2 -2 s 2 


log 




Where £2 is the argument of ascending node and i is the inclination as shown in 
Figure A4. 



Figure A4. Inner Coordinates 
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^ ,AI 2 - 4 “> — (A 12 - 4 J) giv es 

X = SI 


r~ sinhF + A l _ _Q_ , 1 / ?0 - c \ 

n ' ^ + i?? io ® (^j 


ofe 


v _ Ye 


z = 


H sinhF + B 1 _ _Q_ , 1 , /2r>n^ 

25s + i?? los m 


*„wi 


w e __ , 

T* sinh F + C 1 _ Q , _J /20n^ 

2ss+ ^?? iog rf 5 


+ o 



where 


U = U ’ cosft. V Sinn cos r 

V - U' sin 5 + v- cos S cosF 
W = V' sin F 


A =• a ' 


B = A 1 


cos n - B’ sin n cos T 


Sin n + R ' eoo r> t- 

° c °s i 2 cos i 

C = B' sin F ■ 


(A 12 - 45 ) 


(A 12 - 46 ) 


(A 12-47) 


(A 12-48) 

(A 12 - 49 ) 
(A 12 - 50 ) 
(A 12 - 5 1 ) 
(A 12 - 52 ) 
(A 12 - 53 ) 
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U, V and W are the components of V ffl in the reference coordinate system. 


x, e. 


V = Ue. + Ve_ + We, 
— a) —1 —2 —3 


(A 12-54) 


Letting 


L = Ae^ + Be ^ + Ce^ 


(A 12-55) 


and 


R = xe. + ye,, + ze. 
— o —1 —2 —3 


(A 12-56) 


gives 


R. 

— ko 


— sinhF k + t k 
n k 


+ O 


W|s k 


, Q k 1 , | 2Q k n k S k 

" 25 k S k 25 2 S 2 8 \ \ J 


(A 12- 57) 


Replacing sinhF^ by (A 12-23) gives 


1 


*ko (S k> = ^ko S k + ®ko lo « Q k S k + 2 ko + °ko S k lo 8 Q k s 


k ”k 


+ 5 ko S' 1 + F. S' 2 log 2 Q, S. + G. S' 2 log Q. S 


k k 


+ «ko S k 2 + ° ( S k 3 l0 8 3 |S k 


I) 


(A 12- 58) 


where 


-ko -ook 
Q k 

B, = V . 
— ko — cok 


(A 12-59) 
(A 12-60) 
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C 1 
— ko 


Q k, / 2 "k X v XI 
n k \ e k “ rok “ k 


(A 12-61) 


°ko 


3 V -* 

n k 


(A 12-62) 


E, ■ 
— ko 



Q k 

ZlT ^k 

k 


(A 12-63) 


*4co 


Q, 


v 


_3 — cok 
Zn. 
k 


(A 12-64) 


G 


•ko 


Q k l 2 W Q k i 

"k x k/ n k 2n k 


(A 12-65) 


H, 

— ko 



-2 

e k 

+ ■■ ■ V 
+ .—3 -cok 

4n k 


+ i? k log ft) L - 


(A 12-66) 


(A 12-58) gives the behavior of in the overlap domain. The singular 

behavior for large must be matched with the singular terms in the outer 

solution. The matching must also relate the constants and L^. with 

constants of the outer solution. Note that once and are specified the 

hyperbola R^ q is completely determined since n^. and e^ can be obtained from 

V , and L, . 

—cok — k 


A12.2 Second Order 

The exact form of the second order inner solution is given by (A9-7). It is 
an extremely difficult and laborious task to determine the behavior of (A9-7) 
in the overlap domain, i. e. , for large S^. 
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The integral form can be written 


*k2 < S k> 


J B < S k> ^ G k *ko W 


(A 12-67) 


Carlson uses a Taylor series expansion for B(S k , or) but this approach is 
not satisfactory in a second order theory. The largest term in the Taylor 
series expansion is 


B(S k , <r) = I(S k - <r) + o((S k - <r) 3 ) 

From (A12-58) 

§:ko (or) = ^-ko* + 0 ( lo §M> 

C 

Therefore the largest term in is proportional to 
J (S - (r)(rdcr = J (S - <r) - (S - cr)]d (S-<r) 




Since, in general, < S 


< 

/ 


( S - cr) crd<r = + O(S) 


(A 12- 68) 


(A 12-69) 


(A 12-70) 


(A12-71) 
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Therefore the leading term in the expansion of is order and is. easily 

found. The higher order terms cannot be found by this approach, however. 

_ ( 

In the expansion of B the higher order terms are proportional to 

G(R ) (S -V ) n+3 
dS "° 


For n = 0, 1, 2 


dS 


o°<V = G <V 


,1 j dR dR 

^*0’ = = H(Ro> if- 


2 , dR d 2 R 

di 2 + H(Ro) --r 


dS 


dS 


/dR y 


2 

d R 
— o 


dS 


F rom (A 12-58) 


(A 12-72) 


(A 12- 73 ) 


(A 12-74) 


R = A S + 0(log |S|) (A 12-75) 

— o — o 1 

- . -3 

Comparing (A12-72) - (A12-75) with (All-28), (All-36), All-37) and \ 

(A 1 1-89) gives 


C(Ro) = G^S) + 0(log |S|) 


4g(A ) + o(i° gJ S l 
s 3 -° Vs 4 J 


(A 12- 76) 
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H(R ) = H(A S) + O (log |S|) 


lH(A,) 5* - 

-o \ S ; 


(A 12-77) 


T(Ro) = T(AJ3) + O (log |S |) 


T (A ) + O ( lQg J S - 

S 5 -o \ s 6 


(A 12- 78) 


Putting (A 12- 76) - (A 11- 78) into the derivatives of G leads to the conclusion 
that 


,n 


d -G(R ) = O' 1 


dS 


n — o 


n+3 


(A 12-79) 


The higher order terms in B are thus proportional to 
(S - <r) n+3 /s n+3 

Multiplying each of these terms by (A 12-69) and integrating gives the 
contribution of the higher order terms in the Taylor series expansion of 
B, i.e., 


.S JS 

< s - v r - +3 

gn+3 S n+3 

s 


^3 J (S - r) n+3 [S - (S - or)] d (S - o) 


S 2 I 

f s \ 

n+4 _ 

i s z I 

f, V) 

(n + 4) 1 

K 1 ~ s ) 

' (n + 5) 

v 1 ' s) 


n+5 


n + 9n + 20 


+ 0 ( 1 ) 


(A 12-80) 
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Therefore every term in the Taylor series expansion of B contributes at 

2 

least one term of-orde-r S^. Eventually the denominator in (A12-80) would 

become large and after n = N the terms could be ignored. However, in 

2 

order to obtain a reasonable approximation of the contribution to the 
expansion of N would have to be quite large. It will be shown that in 

order to complete the matching R^ must be expanded out to order 
therefore the Taylor series expansion of B is unacceptable. 

What is actually needed is an asymptotic expansion of B^S^, <r) for both 
and c r large. Such an expansion can be shown to differ from the Taylor 
series expansion since no assumption is made that - <r is small. An 
attempt was made to expand ^(S^, tr) to such an order that (A 12-67) would 
be represented by an expansion to order S^. The amout of algebra encoun- 
tered was horrendous and even stymied an attempt to use Formac. It was 
finally determined that such an approach could, at best, only approximate 
the term with an extremely complicated expansion. This approach was 
therefore also considered to be unacceptable. 

Based on (A12-71) and the previously derived behavior of R^ q , (A 12-58), it 
can be assumed that R^ has the following expansion: 

— k2 < S k> = — k2 S k + — k2 S k lo S Q k S k + S k 

+ 5 k2 s k lo ® 2 «k s k + 5k2 s k lo « Q k s k ' 

+ -k2 S k + °( lo B 3 | s k |) ' (A12-81) 

(That the expansion of R^ would have this form was strongly indicated by 
the work done with an asymptotic expansion of BfS^, tr). Some of the terms 
were actually derived but each successive term required an exponentially 
increasing amount of algebra. ) Differentiating (A12-81) twice gives 

2 = 6 — k2 S k + 2 — k2 log Q k S k + 3 — k2 + 2 -k2 

dS, 
k 
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, „ iog Q k S k ‘ 2 D k2 + E k2 ^ (log |S k 

+ 5 « ' S k . S k \ s 2 k , 


Now using (A2-20) and (A 12- 58) gives 


G <SkJ ' = G <* k o S k> + S <*4,0 S k> ^ko lo « Q k S 


■ko k' -ko k “k 


+ « <*ko S k ) G ko + ° « < A -ko S k>- 


1°S S k 
*■ \ K 


* T (A ko S k ) log | S k 


Using (A12-76) - (A12-78) in (A12-83) gives 


G (A ) log Q, S 

G <*W = ~3 S2 ~ + - ( ^ko } -ko 

S k S k 


1 / lo 8 S k 

+ H(A ko )C ko? + 0 (— ^ 

k \ k 


Multiplying (A12-81) by (A12-84) gives 


G(R ko )R k2 = G(A ko )A k2+ C(A ko )B k2 


+ ^ (A. ko ) R. ko A_ k2 s 


!°g Q k S k 


G(A t 


+ H(A ko )C ko A k 2 j_ + 0 «— , 


log 2 |s k 
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Substituting (A12-58) and (A12-85) into (A8-5) gives 


dS 


2 — k2 ~ G lr ^Irn S lr + G \r l0 § G lr S 


k — -ko k k -4c o 


k k 


+ G <^ko> A k2 + G k =ko 


1 log Q k s k 

+ I G( ^o>5k2 + m^ko’Sko ^k2 + G k BfcoJ — 


+ [ G <^ko>5k2' + ffi^ko'^o ^k2 + G k Skol-ST 

L ' J k 


+ Ol 


' l0 8 Tk 


(A12-86) 


Equating functions of in (A12-82) and (A12-86) gives 


^k2 

= G k V/ 6 



(A12-87 ) 

— k2 

= G k W 6 



(A12-88) ' 

— k2 

= [ G <^ 0 >^k2 + G k5ko- 3 ®k 2 ]A 



(A12-89) 

^k2 

= [ G ( A ko >®k2 + H< A ko>S k 2 + G k 5ko]/ 2 



(A12-90) 

*k2 

= GfAfco’SkZ + SAkolSko ^2 + G k 5k„ - 

2 °k2 

J 

IT 

(A12-91) j.rj 


T-oI 

sd" 


Since a PP ears as a linear term in S^ in (A12-81) it vanishes when taking 

the second derivative. Thus it is not possible to get by direct substitu- svs 
tion. However, (A12-81) can be rewritten in the form -oo 

LA) 
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^k2 log Q k S k + -k2 log Q k S k + — k2 


' 5^2 (S k ) / S k " ^k2 S k “ — k2 S k log Q k S k 


5k2 S k + °' 


V| s k | 


(A12-92) 


Letting 


S k = Q k/^k 


(A12-93) 


in (A12-92) gives 


\ 2 log ^ “ -k? l0g + - 


k “ -k2 ^k ' Mc2 


^k Q k ^kZ (Q k/ fi k ) " 


he ^k2 + ^k Q k ^k2 log ^k 


+ ^k Q k ^k2 


+ °(' J k lo 8 3 >*k) 


(A12-94) 


where 


5^2 (Q k/ fX k ) 


/ 


ko 


^k/^k 

B(0 kAV ”■> G k 5ko {r)i ’ 


(A12-95) 


The left-hand side of (A12-94) is exactly the expression which will appear in 

the matching. For very small p's such as in interplanetary applications the 
2 k 

log p^ an< ^ l°g (J-k terms will be much larger than F^, anc * ?4e2 Can Pro^^ly 
be ignored. In such a case (A12-94) is not needed since and can 
evaluated directly from (A12-90) and (A12-91). For the larger value of p^. 
corresponding to the earth-moon system it maybe necessary to evaluate 
(A12-94). This requires first evaluating (A12-95) by numerical quadrature. 
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Then the subtraction indicated on the right hand side of (A12-94) must be 

performed with significant care since it involves the difference of two large 

-2 

numbers, i. e. , the difference of two numbers which are both order . 


The behavior of R^ i n ^he overlap domain is therefore given by (A12-81). 

It has stronger singularities than does but it is also multiplied by pj^ so 

the effect of the singularities is somewhat diminished. The coefficients of 

the expansion (A12-81) can be found as functions of the coefficients in the 

expansion of R, 

— ko 


A12. 3 Third Order 

Only the largest term in the expansion of is of interest. Since the dif- 
ferential equation for R^» (A8-6), is similar to that for it is assumed 

to have a similar expansion. Let - 


*k3 = * k 3 S“ + O(s-Lo g 0 k s k ) 


(A12-96) 


where m is unknown. The second derivative is 


as 2 


= m(m-l)A k 3 S ^- 2 + o(s^- 3 log|S k |) 


(A12-97) 


Multiplying (A12-96) by (A12-84) gives 

G <?W^3 = G '^o>^ 3 s r 3 + °(C 4 l0 «l S kl) 

According to (A8-8) 


(A12-98) 


(A12-99 ) 
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Substituting (A12-98) and (A12-99) into (A8-6) gives 

4> S k + °( S k l0 s| S kl) + °( S r 3 ) (A12-100) 

Comparing (A12-97) and (A12-100) gives 

m '= 4 (A12-101) 

Afc3 = «k ^o/ 24 (A12-102) 

Therefore, in the overlap domain the behavior of R ^ * s g* ven 
Sks'V = h s k + °( S k lo e| S k|) (A12-103) 

A12. 4 Fourth Order 

Like the third order outer solution, the fourth order inner solution is not 
actually used but its behavior in the overlap domain is important. This 
behavior can be deduced without rigorous proof. It has already been shown 
that in the overlap domain 

*ko = °( S k) 1 

*k2 = °( S k) 

= 4k) 

From this sequence it is assumed that in the overlap domain 

Rk 4 = o(s£) (A12-104) 
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A 13 INTERMEDIATE LIMIT 

The overlap domainrhas been defined as the region where t-t k is small in 
the outer solution and is large in the inner solution. In this region a 
new independent variable cr^ is defined by 

°k = ^ t_t pk^/ fJL k’ ° SQ 'o <0 ' <Q 'l ~ 1 (A13-1) 


Note that if a = O, <r, simply shifts the time scale without any scaling. The 

K 

limit a = O is then the outer limit. If n = 1, (A13-1) reduces to (A7-27) giving 
the inner limit. The variable range a <a<a^ will be defined as the intermediate 
limit and a Q and must be defined by the matching. 


From (A7-12) with (3 = 1 



t, + 


P 


k k 


(A13-2) 


Substituting (A13-2) into (A13-1) gives the outer variable 

^~^k = ^k^k ^ ^k T k (A13-3) 

Using (A5-4) 

t-t k = p“ M£<r k + P M k r k (A13-4) 


The last result gives t-t k in terms of cr k and shows explicity that t-t k is 
small, i. e. , order p®. Substituting (A13-1) into (A7-27) gives 


a -1 
^k °-k 


a- 1 

= P 



(A13-5) 


(A13-5) gives S k in terms of o" k and shows explicity that S k is large, i. e. , 
order p ff_ \ since < 2-1 is always negative. 
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(A13-4) and (A13-5) must now be used to transform the outer and inner 
solutions into functions of a' common variable o^. Then thej.match.ing can . 
be carried out. 


A14 MATCHING 

It has been assumed that both the outer solution, (A10-1), and the inner 

solution, (A12-2), are valid in the overlap domain. The matching of these 

two solutions is most simply stated by requiring their difference to be small 

4 

in some appropriate limit. Cole states this limit has 


lim 

p — 0 

cr, constant 
k 


(outer solution) - (inner solution) 


= O 


(A14-1) 


where e(p) is a guage function. For a second order theory c(p) is most easily 
2 

chosen to be p . The limit (A14-1) can be written 


lim 

P-0 

<r, constant 
k 


-ko . -1 


£k 

2 + “ + -2 + ^-3 " ~2 
M- M- 


^k~ko 3 2 4 

+ p Mf R. , + p M, R. - 
p k — kZ k — k 3 


+ ^ M k ^k4 + • • • 


o 


(A14-2) 


In order for this limit to hold all terms which are singular or constant as 

p — 0 must vanish. For the singular terms, the coefficients of the different 

functions of cr^ must vanish identically since the limit must hold independently 

of the value of a, . 

k 


In each of the expansions for £_^ Q , £_p £_y R^, ^k2’ — k3 an< ^ 5^4 

certain terms have been ignored and entered only as order something. These 
terms must all vanish in the limit given by (A14-2). In general these terms 
will vanish in the limit if the exponent of p is positive. That is, p X vanishes 
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if x>0, is constant if x = 0, and is singular if x<0, as ,-*•(). In order to 
make these terms vanish the a of (A13-1) must be restricted as follows 



O 





= O 




5a-2 


Both terms vanish if a> 2/5. 



O 


( 


3 a- 


l ) 


This vanishes if <*>1/3. 





This vanishes if coO. 



The effect of vanishes if «<l/2. 



This vanishes if <*<2/3. 


~ ° 


This vanishes for any a, 


M- 2 ^3 ~ °(h- 2 S l) = Of, 3 - 1 ) 


4a- 1 
F- 
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This vanishes ifa>l/3. 

^ 5k4 ~ °f 3 s k) = °(^ a ' 2 ) 

This vanishes if a> 2/5. All of these restrictions are satisfied if 


? <a 


< * 


(A14-3) 


Therefore 


2 (A 14-4) 

**o 5 


1 . (A 14- 5 ) 

a l ~ 2 , 

It is assumed from this point on that (A14-3) is satisfied and that all of the 
above terms vanish in the limit. Therefore these terms need no longer be 
considered in the limit (A14-2). 

Al 5 INTERMEDIATE FORM OF THE OUTER SOLUTION 


A15. 1 Zeroth Order 
From (A13-4) 


(L-V 

2 

M- 

(tV 

M- 


(tV 

M* 


a - 2 , . , - 1 . , 

H M k <r k + IX M k T k 


2a- 1 .,2a 2 _ a .,1+a , ru,.\ 

» M k °k + 2| " M k T k \ + ° M 


3a-2 ..2a 3 , . Za ~ 1 A/r l+2 a 2 
* M k '"k + 3f " M k T k % 


+ 3^ cr k + O(p) 


(A15-1) 


( A 1 5- 2 ) 


(A15-3) 


89 





4a- 2 ,.40- 4 , So- I 3\ 

I* M k <r k + °l |i V 


+ o( (i 2 “ (T 2 j + o(Ii 1 + ° r k j + o^ 2 ) 


Substituting these expansions into {All- 26) gives 


— ho 2 / _ \ , - 1 

— = A* - Eko) + * 

n ' ' L 


M k T k^ + °^> 


a-2 . .fi 1 , r n .1 + °' * 

"k^ M k M k T k£k 


, 1 a , ,2+0" 2 .. 1+a' 

+ 2 ^ M k T k G k -k + °l tl , 


, 2f 1 2a-l ,,2a * 1 2or-l .,l+2a ' _ 

+ k [" 2 ^ M k ^k + 2 ^ M k k G k 

, ~V 2a\ 3 r 1 3a-2 .-3a ~ .. . „/ 3a- l\ 

+ ;J + °k [z» M k G k + C Y )\ 


. ( 4a - 2 0 4a- 2 -> \ . 

+ 'k V ^ £ V - ^Sr-Sk pj M M k" 


4 4 "- 1 ) 


where the leading term has been carried explicitly in order to show that 
the matching verifies (A10-3) and (All-21). 


A15. 2 First Order 
From (A13-4) 

p _1 log Q k (t-t k ) = p _1 log 2 k + M k ~ a T k ycr k 


i U2 V<) 


(A15-6) 
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(‘-v Q k (‘-v = ^ K \ io 8 z k + ‘4 2 "%, 


+ C>(fi log 2 k ) 


(A15-7) 


1 M k °k + M k \ 


(A15-8) 


-1 ,, , >2 , ~ . . 2a - 1 2 2 . _ 

H (t-t J log Q k (t-t fc ) = [L <T t log s t 


k k 


where 


+ 2|l" <r k log S k 


a , ,1 + 0’ 

+ 11 M k T k °'k 


0(,*-/0 + °(, log = k ) (A15-9) 


S k * ^ Q k M k \ 


(A15-10) 


Substituting these expansions and (A15-2) into (All-78) gives 


£-1 


^ -lk + M k T k -lk + ° (fl) 


+ <r. 


0-1,,°, ' Of ,,0T , 0..1 + 0 T , 

^ M k^lk + ‘ 1 M k -lk + 11 M k k-lk 


, a ,,1 + 0' „ , l+tfll 

2,1 M k T kai k + °V Ij 


+ o-, 


2a- 1 . 2o _ / 2 cm 

^ M k S. lk +°^ jj 
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+ 0 (^)J 

+ log Z fc |> U k c lk + M> + % k Llk + o(^°) 

+ *1 log 2 ^ p"- 1 M^l, k + o(g 2 “) 

, 1 -o 1 -a , _ / 2 - 0 ^ 

+ M k T k^lk + °\'* )} 


+ °(» UZ i<) ... ■; 

(A15-11) 

A15. 3 Second Order 
F rom (A13-4) 


(t-y 1 io g o k (t-t k ) 


= K/\ lo * s k + °(^' Za % z iog \) 

(A15-12) 

<t-, k r' - ^M-»/. k + o(g'-y 2 ) 

(A15-13) 

logQ k (t-y = io g s k+ o(g 1 -'y <rk j 

(A15-14) 


(t-t k ) log 2 Q k <t-y 

-i » ; 

= H* <r k log 2 2 k + O^X log 2 S fc ) 

+ 0 ^ 2 '"/cr k ^ + O^ 2 " 3 * c ^ 2 log ^ (A 15 - 15 ) 


+ log 


-1 

H- ± Y 
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(All-131) can now be written 


£■2 


= [—2k + °<^>] 


+ cr. 


^ M k^2k + ° 


(^°) 


+ log \ 


-2k + 


+ k k log 2 k 


+ k k log Z k 


M k Lk + °(^‘') 


'•*' m£s» + °(«‘ I+ ") 


J k L 


|j. - ^ M7 a b„ + O 
k —2k 




iog 2. 




+ ©(m- 1 Za <r k 2 log 2 fe \ + o([i log 2 2, 


<) + of lo § 2 2 k) 


A15, 4 Third Order 

The effect of r^ is vanishingly small since a<\/2. 

A16 INTERMEDIATE FORM OF THE INNER SOLUTION 


A16. 1 Zeroth Order 
F rom (A13-5) 

log Q k S k = l0g 2 k " l0g ^k 


(A15-16) 


(A16-1) 
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(A12-58) can be written 


M R, y l ■ 

= >*• ' K ^ko \ + M k ?ko < lo S A - ^ ^k> 


1 Q- ?.« ( lQ g 2 k ‘ l0g fX k ) 

+ fi" 1 M. C, + hl M“ D. — 

r k — ko n k — ko o r 


+ K' a ?ko k + 0 (^' 2 “ "k 2 l0 8 2 \) ' < A16 - 2 > 


Al6. 2 Second Order 


(A12-81) can be written 


, ,3 3»-2 . a 3q . 3 

IX M k 5^2 = M k <r k 


2n-l . rl+2n „ . 2 

+ M- M w B (log S t - log n k ) cr k 


, 2ar- 1 , ,1 + 2(2 „ 2 

+ •* “k C -k! \ 


+ M k + ° 5^2 ( l0 B 2 2 k -2 log |x k log E k 

+ lo e 2 l‘k) <r k + ^ M k + “ ?k2 <lo 8 2 k 


1 \ 1 ® 71^2+°' T' 

- lo § *V °k + >\ M k ^2 °k 


(A16-3) 


A16. 3 Third Order 
(A12-9&) can be written 


fx 2 R. 


k — k3 


4a- 2 
M- 


•» , 4 i/ A 

M k ^k3 


(A16-4) 
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A16. 4 Fourth Order 


The effect of is vanishingly small since o< 2/5. 


A17 RESULTS OF THE MATCHING 

From each of the terms in (A14-2) similar functions of can be collected 
from (A15-5), (A15-11), (A15-16) (A16-2), (A16-3) and (A16-4). Since the 
limit in (A14-2) must be satisfied independently of the value of the 
coefficient of each function of cr^ must vanish. The following results are 
then obtained. 

(a) Terms proportional to <r^: 


o-2 


M- 




o 


M 


l + o 


k T k Ek + 2 11 


o .,2+o 2 ^ 

M k T k G k^k 


o-l ..o , , o . , o , a . , 1 + o , 

+ p M k d lk + p M k e lk + p M k T k f_ lk 


, o .,l + o Q a ( l + o\ 

+ 2( * M k \£,k + ' 1 M k&2k +C T ) 

o-2 . .a . o . ,2+o 2 

M k ^ko - ^ M k 5k2 lo 8 ^ k 


o 2+0 
-p M E 


-4c 2 


i o a / t 2+ q ' TT 

log^k-l" M k ^2 


= o 


(A17-1) 


or 


+ 1* 2 m2 ( a ,? io? 2 - a .2 lo « ^ + ?*2 


\ + ^Ik + F (£ lk - ™ k \ Ek + \ M k i G k ¥k 


+ M k + 2M k T k i-lk + &2k + W > 


(A17-2) 


(b) Terms proportional to log 2 • 

1C 

^ _1 a 1]c + c lk + O(p) - p- 1 M k B ko = O 


(A17-3) 
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or 


5™ = + £7> + °(t‘ 2 » 


(A17-4) 


(c) Constant terms: 


l*' 2 (iko - Eko) + 1 M k V k E |T 1 b lk + M k T k d lk 


+ i 2k + °M * l *" 1 M k Ifeo lo § H k • I *" 1 M k ^ko = ° < A17 - 5 > 


or 


^ko " — ko ^ ^-^ok " ^-ko) 


+ \ Yfc + M k X kjk 


+ ^( T k -lk + M k -2k) 


+ 0(p 2 ) 


(A17-6) 


It is clear that in order to balance the left and right hand sides of 
(A17-6) that the order p - * term must vanish, i. e. , 


^ok ^ko 


= O 


(A17-7) 


(A17-7) is equivalent to (A10-3) and (All-22) and verifies the 
earlier assumption. Then (A1 7- 6) become s . 

5kO = ?ko lo « “k + T k + M k‘ tik 


+ ^( T kilk + M i I ^k) +0( ' l2) 


(A17-8) 


(d) Terms proportional to cr^ 1 log 2^: 


-a . .-o , rs , -os -a , .2-0 

M k —2k + 0(|X ) ‘ ^ M k 5ko = ° 


(A 17-9) 
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or 


2ko = M k 2.Zk + ° (|j) 


(A17-10) 


(e) Terms proportional to 


-a , ,1-Q' , -a . ,-a , , r 1-a. 

>* M k T k 2-lk + 11 M k —2k + 0(11 > 


-O' „ r Z -a 


-<* w 2 -® 


+ * M k ^ko lo S ^ £ 


k — ko 


O 


(A17-11) 


or 


5ko = lo B l*k + M k' T k 2.1k + M k 2 —2k + ° M 


(A17-12) 


(£) Terms proportional to o^: 


1 3a-2,,3a_ 3o- 1 . 3a-2-..3a. — 

M k G k ^ + °(p ) - |x M k A^ z = O (A17-13) 


or 


k2 = Z °k — k 


G, V., + O(p) 


(A17-14) 


(g) Terms proportional to cr log E • 


2o-_l 2 a . 2a. 2a-l l + 2a „ _ 

M- M k _f lk + 0(p ) - p = O 


k t -k2 


(A17-5) 


or 


^kZ “ iVi k -lk 


+ O(n) 


(A17-16) 
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(h) Terms proportional to cr£: 


1 2a- 1 *„2 a 


" 2 ^ 


M 


k £k + 2^ 


1 2a- 1 


1 + 2 * 


M, 


\ G k ^k 


or 


2a- 1 .,2a ~ 2a. 2a- 1 . .l + 2a„ , 

+ H M k g. lk + 0(|j. ) + p M k B k log p k 


2a- 1 ^,l+2a ~ ~ 

-p M k = O 


^k2 " — k2 log ^k “ 2 M k £k + 2 \ G k ^4c 


r~l 


+ M k £ lk + OM 


(i) Terms proportional to <r k log 2 k : 


o ,-.,1 + 0. a ,,2+ a - 

M- M k £ 2k + 0(p ) - p M k D^ z = O 


or 

^k2 = M k -2k + ° ( ^ ) 

(j) Terms proportional to <r k log 2 k : 

a , o a .,1 + a r , o u tt c 

V- M k £ik + 2fX M k T k— lk + ^ M k -2k 

. l + o. a .2 + a -t- . 

+ 0(p ) + 2p M k lo S H-k 

a . ,2+o - 

' 11 M k ?k2 = ° 


(A17-17) 


(A17-18) 


(A17-19) 


(A17-20) 


(A17-21) 
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or 


" 2 l0g ^ + —it + 2 M t T t— 1 


^kZ -5 r k • — k _i k 


k k— lk 


+ M k 2 l2k + ° ( ^ ‘ 


(k) Terms proportional to cr^j 


1 | 4a-2 -,4a 2 . > 4a-2 . .4a u • 2 . . 


24 


. 4a-l. 4a-2 ,,4a 

+ 0(p ) - p M k = O 


or 


^k3 


Sk 


^ (t k> - (t k> 


h 


24 + O(p) 


(A17-22) 


(A17-23) 


(A17-24) 


(1) The terms not accounted for so far are of the following orders: 

4a- 1 l-2a 
P , p , P 

These terms vanish identically in the limit p— -0. 

A18 SOLUTION OF THE INITIAL VALUE PROBLEM 

The expressions presented in Section A18 contain the solution to the following 
problem: given a set of initial conditions at t = t what are the parameters 
defining a close approach of the kth body at or near t = t^? 

r(t 0 ) = £ 0 (t 0 ) + P£ 1 (t o ) + p 2 £ 2 (t o ) (A6-13) 

v(t Q ) = v Q (t o ) + p Vj(t o ) + p 2 v 2 (t Q ) (A6-14) 

Given these constants plus t and t, defines all the constants of the outer 

r o k 

solution (including integral constants). 
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The constants of the inner solution are the vectors R_ (S. ) and V, (S, ) 

—ko ko' — ko v ko' 

which appear in (A9-1). The higher order terms are set to zero in (A9-6) 
meaning that at = S^ 0 the zeroth order inner solution (an hyperbola) and 
the second order inner solution (a perturbed hyperbola) are equivalent. 
Therefore the complete inner solution is determined at if the corre- 
sponding zeroth order hyperbola is known. 


Rather than determining (S^ q ) and (S^ q ) as functions of rjt Q ) and 

v(t ) it is more convenient to determine the elements of the zeroth order 
— o 

hyperbola. They are 


a 


e 




semi-major axis 

eccentricity 

inclination 

argument of the ascending node 
argument of pericenter 
time of pericenter passage. 


(A18-1) 


Usually it will be desired to determine the inner solution at closest approach. 
Then = 0 and the elements of the zeroth order hyperbola can be used 
directly to calculate the closest approach distance, etc. .. . . 


In order to determine the relationship between the set of elements (A18-1) 

and the initial position and velocity it is necessary to rewrite (A17-2) and 

5 

(A17-8) in a new form. From Battin the unit normal to the plane of motion is 


*k 


= ( sin sin i, , 

k k 


- cos n sin i^, 


cos 1, 


(A18-2) 


The impact parameter vector is defined as the vector normal to Y-cok with 
magnitude equal to the semi-minor axis b where 


b 


k 



(A18-3) 
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Figure A5. Impact Parameter Vector 


For the trajectory shown in Figure A5 the unit normal N^_ is out of the plane 
such that b^. and form a right handed system.' Then 


V 


K 


= b 


c°k 


k V, 


x N. 
— k 


(A18-4) 


Substituting (A12-54) and (A18-2) into (A18-4) and then using (A12-38) - 
(A12-39), (A12-48) - (A12-53) and (A12-55) eventually gives 


Q k 

b, = L, + =— V m , . 

— k — k n. — ®k 

k 

i 


(A18-5) 


where 



Then (A12-61) becomes 



- 1 


V m , + b, 
— ®k — k 


(A18-6) 


(A18-7) 
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V 


The expressions for A ^ p and C^, (Al'7-2) and (A17-8) become 
-«=k = A + ^-lk + 11 1 k + &2k + M k T k (-lk + 2 &lk 

- Ek) + I 

- ^k2 lo « l"k + ^ 2 )] + ° (ti3) 


M k T k G k ^kf M k(5k2 lo s 2 •*. 


Q, 

-k " T k -«k ~ h. 


2n, 


log 


^ 1 , e 


-1 


k k 


V 


“k 


(A18-8) 


+ M' 1 (b lk + pd_2 k > + 0(p 2 ) 


(A18-9) 


Since (A18-9) is only accurate to order p. the p term in can 

temporarily be ignored. Also, since b^ is normal to V , taking the l nner 
product of (A18-9) with gives 


Q, 


n, 


2n. 


lo g| 


^k e k; 


-1 


JC (b lk + ^ 2k ) • v„ k + 0(^ 2 ) (A18-10) 


(A18-11) 

where I is the unit diagonal matrix and the transpose of V^. The 

transpose is introduced through the identity 

(x : y)z = (l l T )x (A18-12) 

, , ; \\ { \''\ ■ \ \ 

The solution is now possible.. The^ first two 1 tefms (of (i^il8-8) are used td ■ 

evaluate V . . Then 1 

— ®k 


Subtracting this result from (A18-9) gives 

' T \ 

M. \‘ ~ _2 / ( -lk + ^-2k 


*k 


1 





(A18-13) 
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and n^_ is found from (A18-6).. .Next (A18-11) is used to determine and 

\ - N ,AI8 - 14) 

®k - ( 1+ ^Ak)’ /2 ' " (A18-15) 


Since £^ k and i^ are as yet unknown, the unit normal is determined from 


Nk = ^ k *v <ok )/(B k v IBk ). • . . 

(A18- 1 6 ) 

Then from (A18-2) 

. t 

cos i^ = • e^, OS i^ - w 

(A18-17) 

sin = (N k • e^/sin T k 

(A18-18) 

cos ^ k = - (N k • sin i 

i 

(A18-19) 

Finally, inverting (A12-48) and (A12-49) gives 

<J k = U k cos ^ k + V k sin <T^ k 

(A18-20) 


V k = ( v k cos f2 fc - U k cos ^/cos i (A18-21) 




(A18-23) 
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Enough information is now available to evaluate (A12-94) and this, along 

2 

with (A18-10), allows the order p. term in (A18-8) to be evaluated. 


The time of pericenter passage is found from (A13-2), i. e. , 


* P k = ‘k + •* T k 


(A18-24) 


The complete set of orbital elements defining the zeroth order hyperbola 
can be evaluated from the expressions for A ^ and obtained through the 

matching. The remaining expressions obtained from the matching in 
Section A17 give redundant information which can be used as an analytical 
check on the matching process. For example, from (A12-60) 


-ko ~ " Q k- ( -®k^ 


(A18-25) 


After some manipulation, (A17-4) can be reduced to 


B ~ - Q, 
— ko k 


KY^ + p G(V±) d lk + 0(/) 


(A18-26) 


Equating (A18-25) and (A18-26) gives 


l(v« k ) = KYk) + > d lk + 0(p ) 


(A18-27) 


Comparing (A18-27) with (A2-13) and (A2-14) leads to 


V„ k = V* +> d lk + 0 (l x‘) 


(A18-28) 


Since (A18-28) is contained in (A18-8) no new information has been, obtained. 


As another example, (A12-62) gives 


D - V V 
~ko -«k/ V °°k 


(A18-29) 
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while (A17-10) can be reduced to 


Ac = i/ V k +0 ^ 

using the identity 

G < V A * 2 W v k 

Comparing (A18-29) and (A18-30) gives 


(A18-30) 


(A18-31) 


Y^ k = + O(p) (A18-32) 

a result which is also contained in (A18-8). 

It can also be shown that reproduces information contained in C^. The 

proof involves replacing by b^ using (A18-5) in (A12-61) and (A12-63) 
then taking scalar products with and b^. The algebra is long and tedious 

and requires certain identities such as 

G(x) y x =-^-(x. y) (A18-33) 

x 

Since the algebra is so lengthy and complex no further proofs will be given 
here but only the statement that such proofs have been worked out and the 
matching appears to give a consistent set of equations relating the constants 
of motion of the inner and outer solutions. 


Although the notation is different the results obtained by Carlson are contained 
in (A18-8) - (A18-11) except that he has one order less accuracy in each 
expression. In terms of a common notation the results are identical to first 
order for and to zeroth order for and b^. The idea of decomposing 

(A18-9) into orthogonal components to solve for and b^ independently 
comes from Carlson who introduced b^ directly into the development of the 
inner solution prior to matching. In Section A12 the vector L, was introduced 
and, from (A18-5), 


^k 


4 


Q k 

IrT^k 


(A18-34) 
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The representation of in terms of b_^ and is certainly advantageous 
for solving the initial value problem, as presented in this section, but it 
is not unique. Other representations may also be useful, in particular for 
two-point boundary value problems. Therefore the solution developed in 
this study and contained in the expressions for and in the preceding 
section is somewhat more general than that of Carlson as well as being of a 
higher order. The development of the boundary value solution from and 
appears in Section B. 
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1 > Section B 

SECOND ORDER TWO -POINT BOUNDARY VALUE SOLUTIONS 

B 1 FUNDAMENTAL SOLUTION 

The fundamental relationships between the constants of the outer and inner 
solutions in the second order asymptotic solution of the problem of N bodies 
were derived in Section A17. In Section A18 it was shown how these equations 
can be used to formulate an initial value solution. In this section it will be 
shown how these same equations are used to formulate certain boundary value 
solutions of practical interest. 

The fundamental equations resulting from the matching were shown in Section 
A17 to be (A17-2) and (A17-8), They can be rewritten as 
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; whe re 

<: ^k :.iik- A <vyiA>- B ivyiiiy < Bi - 3 » 

-Ik = -lk ' C<t k’ ‘o' -l^o* ' D(t k' t o ) — l (t o’ (Bi - 4 > 

4k = —2k - A < t k t o»i2< t o> - B ' t k-'o»i2< t o> ' B1 - 5 ) 

■£ k .=. - C <V l o> %<‘o» - *><***.» ■' (B1 - 6 » 

^2 = '' M k(-Ek2 log " — k2 log g k + — k2 ) (Bl-7) 


Given initial conditions along the outer solution, (Bl-1) and Bl-2) give the 
values of , C ^ and t^, the constants of the inner solution. However, 
using (A3-13) and (A3-31), (Bl-1) and (Bl-2) maybe inverted to give the six 
component state vector. 




(Bl-10) 





(B 1 - 1 1 ) 
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% 


* 

'&2k 


" -lk 


— k2 


“ M k T k 


(£ik + 2 & 


lk ^k 


2 M k T k G k- 


i) 


(B 1 - 12) 


Given initial conditions along the inner solution (Bl-8) gives the perturbations 
to the zeroth order ellipse at t = t . Thus (Bl-8) is an inverted form of (Bl-1) 
and (Bl-2). Either form may be considered as the fundamental solution from 
the matching. 

B2 ASYMPTOTIC BOUNDARY VALUE SOLUTIONS 

The fundamental solution from the matching is essentially an initial value 
solution. However, with a little algebraic manipulation the solution can be 
transformed into a boundary value solution where some variables are pre- 
scribed at each end of the trajectory and other variables at each end are 
dependent on the prescribed values. The goal then is to write (Bl-8) in such 
a form that the dependent variables appear as explicit functions of the inde- 
pendent variables with the independent variables chosen to represent a 
realistic boundary value problem. 


It is first necessary to derive some expressions for the inner hyperbola. 

Suppose the excess velocity, _V t ^ le inclination, i^> the pericenter radius, 

p, , and the time of pericenter passage, t . , are assumed to be known for a 
K th pk 


close approach to the k 
defined by 


body. The excess velocity has cartesian components 


^cok 


= (U 


k' v k’ w k> 


(B2-1) 


V . in 
— ook 

the (x, y, z) 

re B 1 


l>* 

II 


II 

y.^ir 

1 ^ 

II 



(B2-2) 

(B2-3) 

(B2-4) 
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Figure Bl. Hyperbolic Excess Velocity 

Since the plane of motion must contain V_ t ^ le ascending node must satisfy 
(cf Battin 5 , P. 179) 

?2 k = \ k ± F k + (1 ± 1) tt’/Z (B2 


where 

r • 

(T, = tan O', /tan i. , 0 S or < -n- /2 (B2 

k k k k 

With some manipulation (B2-5) becomes 


/ — 2 — 2 

-k - \ u k + v k 


-1 


cos n, = ur + vr J ctn i, \ V, W. tU, u: + vr tan i. - w 


1 


■kl-k-'k" “kl( u k + v k) u p 2I k- w kf 


(B2- 


■Yu? + v? ) ctn 


iiju, W, ± V, f(u? + V?' 


tan 2 i -.W 2 | 2 1 (B2 - 


H 


sin £2 



It is obvious that 


ta ” 2 r k 2 w k / < u k + v k» 


(B2-9) 


giving a minimum i^. which is compatible with (B2-5), (B2-7) and 

(B2-8) the upper sign is to be used if the approach or departure is to be over 
the kth body (b k • e_ > 0) and the lower sign is to be used if the approach or 
departure is under the body (b^ • e_^ < 0)* 

The pericenter radius is defined as 


p k = a k (e k - *> 


where the semi-major axis is given by 


a k V ook 


(B2-10) 


(B2-11) 


Then the eccentricity is 

p k = 1 + p k v .k 2 < B2 - 12 > 

From (A18-22) and (A18-23) 

sin = -r k I/2 [(5 2 - l) 1/2 u; + Q k v;] /e k (B2-13) 

cos = a' /2 [(? 2 - l ) 1 12 V ' - Q k u ' )/ e k (B2- 14) 

where, from (A18-20) and (A18-21) 

Uk = U k cos + V k sin 

\ \ * 

V k = ( V k cos n k " U k ® in / cos i k 

= ^k /sin ^k (B2-17) 


(B2-15) 

(B2-16) 
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Equations (B2-7) - (B2-17) define the orbital elements (in addition to i^ which 
is assumed known) for the inner hyperbola. From (A 12 -27) and (A 12 -28) 


A k = a k e k cosw k 


(B2-18) 


B k = Vk Sinw 


C'-ii f. 


(B2.-19) 


and from (A12-51) - (A12-53) 


Ar = A^ cos - B^ sin cos i^ 


B, 


A^ sin + B^ cos cos i^ 


B k sin i k 


(B2-20) 

(B2-21) 

(B2-22) 


Then, from (A12-55) 


= (A k’ B k' C k> 


(B2-23) 


The two vector constants of the inner hyperbola are given by (A12-59) and 
(A 12 -6 1 ), i.e., 


^ko 




c i 

— ko 


Q k / 9 nk^ 

L +— -log(^iS) V 

~ k n k & \ e k/- 


®k 


whe re 


n, 


= V 


cok 


(B2-24) 


(B2-25) 


(B2-26) 


The equations of this section show that known values of *k an9 ^k are 

sufficient to define the inner hyperbola and its two constants, A^ q and C_ ko * 
It is now possible to proceed with the boundary value solution. 
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B2. 1 Midpoint-to-Target Body Solution 

The basic midpoint-to-target body solution is shown in Figure B2. The initial 

time, t , and position, r(t ), as well as the pericenter time, t , radius, p^, 

. o — o pi 1 

and inclination, i^, for a close approach to the target body k = T are all spec- 
ified, The initial velocity at t = t is unknown and must be determined. An 
ephemeris is required giving positions and velocities of the N-2 perturbing 
bodies with respect to the primary body. The ephemeris fixes the coordinate 
system of the asymptotic solution. 

The initial conditions for the asymptotic solution are given by (A6-13) and 
(A6-14). They are 


£ (t o } = ^o (t o } + ^l (t o } + 11 


(B2-27) 


^ (t o } = ^o (t o ) + + H- — 2^ t o ) 


(B2-28) 


From Figure B2 


r (t ) = r(t ) 
— o o — o 


(B2-29) 


the refore 


^l (t o ) = — 2 (t o ) = ° 


(B2-30) 


For simplicity let the initial velocity perturbation be defined by 6v(t Q ), i. e. , 


6 ^(t 0 ) = li(t 0 )+ ^ 2 (t 0 ) 


(B2-31) 


From Figure B2 the final position of r^t) is 


~ ET (t T ) 


(B2-32) 


113 




where, from (A13-2) 


fc T t pT " t T 


(B2-33) 


and is the position of the target body obtained from the ephemeris. The 
parameter t^, is not prescribed by the boundary conditions and can, without 
loss of generality, be set equal to zero. 

The two position vectors given by (B2-29) and (B2-32) define a Lambert 

problem and the solution gives r^(t), shown as the dashed line in Figure B2, 

and the velocities, v (t ) and v (t_). 

— o o — o T 

Now let k = T and Q^, = -1 (cf. (All -34)) with 


m t = 1 


(B2-34) 

(B2-35) 


and substitute (B2-24), (B2-25), (B2-30) and (B2-31) into (Bl-8). Using 
(A3-3 1) gives 


B(t T ,t Q ) 6v(t o ) = L t - 


1 


P T e T 


T T ' ”t 108 V^T 


- 1 , 


E)(t T , t Q ) 6v-(t Q ) = p (V a)T - V T ) + 6 T + pjl T 
Solving for bv(t Q ) and gi yes 


X^t + X T + H-£ T (B2-36) 


(B2-37) 


6v(t o ) = B(t T ,t o ) 


-1 


+ V_rp + 


— T 

1 . /^t 6 t\ 

r T " TT t log \ 2n r _ j 


T \ T / 


— ®T 


— 00 T — T + ^ D ^T’ t o' > — ^o^ ” T ^ 2 -T 


(B2-38) 

(B2-3^) 
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These two expressions are not explicit since V appears in (B2-38) and 
6_v(t ) appears in (B2-39). However, if (B2-39) is substituted into (B2-38) 
an expression involving only 6v(t Q ) results. Rather than solving such an 
expression for 6v(t Q ) the solution can be obtained as follows: 

1. Solve the Lambert problem defined by (B2-29) and (B2-32).. This 
gives ^ 0 (t Q ) and _v_ q ( t T ). 

2. From (All-22) 

V T = v o (t T ) - £ T (t T ) (B2-40) 

where £^,(t^,) is obtained from the ephemeris. 

» ' - • 

3. Let V _ = V^,. This gives the zeroth order excess velocity. 

4. Evaluate (B2-1) - (B2-26). 

5. Evaluate (B2-38) using the results of step 3 and 4 and with 4^ = 0. 

This gives the first order velocity correction. 

6. Evaluate (B2-39) using the results of step 5 and with r^, = 0. This 

gives the first order excess velocity. , i 

7. Repeat step 4 using the first order excess velocity. 

8. Evaluate (B2-38) using the results of steps 6 and 7. This gives the 
second order velocity correction. 

9. Evaluate (B2-39) using the results of step 8. This gives" the second 1 
order excess velocity. 

10. Using the results of Steps 1 and 8 in (B2-28) gives the second order 
initial velocity 

v.(t 0 ) = _v 0 (t 0 ) + li&vjto) (B2-41) 
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11. Repeat step 4. This gives second order orbital elements at the 
target body. >■ 

Throughout all the steps the prescribed values of r_(t Q ), P^,, i^ and t ^ remain 
constant. Steps 10 and 11 represent the second order solution to the boundary 
value problem. 

B2 . 2 Launch Body-to-Target Body Solution 

This solution is shown in Figure B3. The pericenter time, t radius, p^, 
and inclination, i^ are prescribed at both the launch planet, k = L, and the 
target planet, k = T. The hyperbolic excess velocities at both planets are 
unknown and must be determined. An ephemeris is required giving positions 
and velocities of the N-2 perturbing bodies with respect to the primary body. 
The ephemeris fixes the coordinate system of the asymptotic solution. 


From Figure B3 it can be seen that the zeroth order ellipse passes through 
the launch body at t = t and through the target body at t = t„, where 

i-i 1 


= VVk ,B2 ' 42) 

where again can be set equal to zero without loss of generality. Then 

£o(t L ) = £ L (t L ) ' (B2-43) 

r o (t T ) = £ T (t T ) (B2-44) 


where £ and £^, are the positions of bodies k = L and k = T obtained from the 
ephemeris. The two position vectors given by (B2-43) and (B2-44) define a 
Lambert problem and the solution gives r^(t), shown as the dashed line in 
Figure B3,: and the velocities, v (t T ) and v (t„). 

Now let k = L, T with 


P = or p T (B2-45) 

Q l = +1 (B2-46) 
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Target 

Body 
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Figure B3. Launch Body- to- Target Body Solution 


-1 


(B2-47) 


^cf (A 11 -34)^ and substitute (B2-24) and.(B2-2 5) into (Bl-8) for each k. This 
gives two equations each of which have the same left hand sides. Assuming 
the position and velocity continuous at t = t allows the two right hand sides to 
be equated ^ Using (A3 -3 1) and (A3 -32) gives 


A ( t T’ 'l.) M I>L- M L T L + ^ lo s(T^) — coL + -L + ^-L +B < t T- t L ) . 


^ coL " — L* + -L + 


i , r 1 T e T 


= M t L t -M t t t -^-logl^l 


+ Y_rj> + |lC_rp 


(B2-48) 


c(t T , t L ) m l l l - m l L ^iog(^) v mL + Y L + p£ l + D(t T . t L ) 


• - X L ) +1 L + - X T ) + 6 t + ka T (B2-49) 


Solving for V L and V^ T gives 


— cdL _ — L " ^-L “ Hr + H- B ( fc T’ ^ M t— T ' M T T T 


- J- log ( -£-£)! v_„ + + fi£ T l -pB(t T ,t L )" A'(t T .t L# 


2h T ]\ — ®T ' -T n -T 


i M I>L- M L | T L + h 7 l0g coL + -L + h -L i 
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X ,«p “ H~ “ H- n T 


M L^ L - m l 


t t + log 

L S L 


>L e L' 


2n 


L/J 


-coL + — L + ^ 


+ D(t T , t L ) 


-coL - V L + I^A L + 1^ H L 

(B2-51) 


If the excess velocities are used to calculate orbital elements then the position 
and velocity at pericenter are defined by the well known expressions 


R , = ^r^r- L. 

— pk a k ^ k — k 


(B2-52) 


1 I 1 + e k\ 1 /2 8 — k 


-pk 


a k 6 k \ Pk 


3 co . 


(B2-53) 


The two expressions for the excess velocities are not explicit but the solution 
can be obtained as follows: 

1. Solve the Lambert problem defined by (B2-43) and (B2-44). This , 
gives (t j_, ) and 

2. F rom (All -22) with k = L, T 

L = A, V - i k (‘ k ) < B2 - 54 > 

where ft. (t^ ) is obtained from the ephemeris. 

3. Let V . = V. , k = L, M. This gives the zeroth order excess 

—CO K K. - 

velocities. 

4. Evaluate (B2-1) - (B2 -26) with k = L, M. 

5. Evaluate (B2-50) using the results of steps 3 and 4 and with 

= tj T = 0. This gives the first order V T . 

— L — co L 
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6. Evaluate (B2-51) using the results of steps 3, 4 and 5 and with 

- t, T - n T =0. This gives the first order V "•••' 

— — Lj — CO I 

7. Repeat step 4 using first order excess velocities. 

8. Evaluate (B2-50) using the results of steps 5, 6 and 7. This gives 
the second order V T . 

9. Evaluate (B2-51) using the results of steps 7 and 8. This gives the • 
second order V. 

10. Repeat step 4 using second order excess velocities. This gives 
second order orbital elements. 

11. Evaluate (B2-52) and (B2-53) using the .results of step TO. This gives 
second order positions and velocities at the pericenters of the hyper- 
bolic trajectories close to the launch and target bodies. 


Throughout all the steps the prescribed values of t and i^, k =' L, M 

remain constant. Steps 10 and 11 represent the second order solution to the 
boundary value problem. 


Although the midcourse time t Q does not appear explicitly in (B2-50) and 

(B2-51) it does enter implicitly since all of the constants are evaluated 

either between t T and t or between t and t_. The midcourse time is not 
-Li o o 1 

fixed and may be defined by 


fc o = {t L + t T )/2 (B2-55) 

The first order position and velocity perturbations at t = t are needed to 
evaluate the second order constants. They are obtained directly from (Bl-8), 
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B2.3 Non-linear Solutions 

Evaluating (B2-38) and (B2-50) requires the inversion of the B-matrix, This 
inversion may tend to give inaccurate results in cases where the linear state 
transition matrix is not a good approximation due to large non-linear effects. 
Such cases arise if one of the endpoints is close to pericenter or apocenter of 
the zeroth order ellipse. The form of (B2-36) (B2-37) (B2-48) and (B2-49) 
suggests an alternative approach. This approach has been called the non- 
linear solution by Carlson. 3 


The zeroth order ellipse is defined by its position and velocity, r^t) and v^(t). 
Suppose a neighboring ellipse is defined by 


r^(t) = r^t) + pAr^t) 


(B2-57) 


v^(t) = ZqW + 


(B2-58) 


where Ar (t) and A v (t) will be called the offset position and velocity. Since 
— o — o 

r^(t) and v^ (t) also define a two-body trajectory the offset positions and 
velocities at any two times t^ and t ^ are related by 






= •< t 2' t l ) 






(B2-59) 


122 



Now define and Y_ k by 


*k = - M k 


Q 


k , /^kV 


T k + H” log \ 2- 


nk /J 


¥,k + \ + ^ 


( B ?. - 6 0 ) 


Xk - ^ (I fflk - Vfc) + Ik + ^-k 


(B2-61) 


Then the midcourse solution, (B2-36) and (B2-37), becomes 



4>(t T . t Q ) 



(B2-62) 


while (B2-48) and (B2-49) become 



(B2-63) 


Now note the similarity between (B2-59) and (B2-62) and (B2-63), This sug- 
gests that (B2-62) and (B2-63) represent the propagation of offset end condi- 
tions for a new zeroth order solution. 


For the midcourse solution, comparison of (B2-59) with (B2-62) gives 


A Io (t o' = 0 

(B2-64) 

A £ 0 ( t T ^ = —T 

(B2-65) 

Then from (B2-57) 


r' (t ) = r (t ) 
— o o — o o 

(B2-66) 
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The end points r^(t Q ) and r^(t^,) define a new Lambert problem, the solution 
of which gives r^(t), shown as the dashed line in Figure B4, v_^(t Q ) and v^(t ). 
Again comparing (B2-59) with (B2-62) gives 


Av (t ) = 
— o' o 

6 ^ t o ) 

(B2-68) 

A -X> (t T ) = 

X T 

(B2-69) 


By combining (B2-41), (B2-58) and (B2-68) the initial velocity becomes 

v (t ) = V' (t ) (B2-70) 

— o — o o 

By combining (B2-58), (B2-61) and (B2-69) the excess velocity at the target 
body becomes 

■ X„T = iy- K-6 T -Va x (B2-71) 

where 

Xt = X (t T ) ' £ T (t T } (B2-72) 

Equations (B2-70) and (B2-71) represent the non-linear solution. They replace 
(B2-38), (B2-39) and (B2-41) of the solution discussed in Section B2. 1. 

For the launch body-to-target body solution comparison of (B2-59) with (B2-63) 




Then from (B2-57) 


L'V = W + t^L (B2-75) 

r^(t T ) = £ c (t T ) + nX T (B2-76) 


The end points r' (t ) and r 1 (t ) define a new Lambert problem, the solution 
o JL/ o _L 

of which gives £^(t), shown as the dashed line Figure B5, jv^(t ) and v^(t^). 
Again comparing (B2-59) with (B2-63) gives 





(B2-77) 


^V 0 (t T ) = Y t 


(B2-78) 


By combining (B2-58), (B2-6 1), (B2-77) and (B2-7 8) the excess velocities 
become 


^L 

= IL - 

(B2-79) 

— ooT 

2 

= Y-T ' “ k 1 T 

(B2-80) 

where 



II 

- e-M 

(B2-81) 

n 

v^(t T ) - E T (l T ) 

(B2-82) 


Equations (B2-79) and (B2-80) represent the non-linear solution. They replace 
(B2-50) and (B2-51) of the solution discussed in Section B2.2. 


Since X. is a function of V , the non-linear solutions, like those of Sections 
— k — cok 

B2. 1 and B2. 2, must be evaluated in a sequence which uses the best available 
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th 

approximation for V , i.e., the n order solution for V . requires the 

fcl"l COxC coK. 

(n-1) order value of to evaluate the offset endpoint X^. 

The non-linear solutions require more computation time than the linear or 
standard solutions since solving a second Lambert problem is, in general, 

i 

more time consuming than inversion of the 3x3 B -matrix. 


B3 APPLICATIONS OF THE BOUNDARY VALUE SOLUTIONS 

The solutions of Section 3 can be used to solve several boundary value prob- 

i 

lems. They are discussed in the following sections. 


B3. 1 Earth-to-Moon 


The simplest earth-to-moon boundary value problem is shown in Figure B6. 
The initial time, t , the initial position relative to the earth, rjt^), and the 
pericynthion radius, p^, inclination i^, and time, t ,p, are all prescribed. 
The initial v 
determined. 


The initial velocity relative to the earth, v(t Q ), is unknown and must be 


This problem is solved using the midpoint-to-target body solution of Section 
B2. 1. The primary body is the earth and an ephemeris is required giving 
the motions of the moon, sun and any other significant bodies in cartesian 
coordinates with the earth at the origin. Although the sun's mass is large 
compared to the moon's mass, its effect is diminished by its great distance 
from the earth. As discussed in Section A1.2 both the moon and sun contrib- 
ute effects of order p where p is the dimensionless mass of the moon. 

In this problem the subscript T of Section B2. 1 refers to the target body which 
is the moon. The effects of the sun and any other bodies enter only through 
the constants V_ T> 6_ T , and rj_ T . 

B3.2 Earth-to-Moon Midcourse 

In the previous section the initial position, r_(t o ), was implicitly assumed to 
be close to the earth. The same analysis may also be used for a midcourse 
maneuver where the position, r(t ); represents a point between the earth and 
the moon as shown in Figure B7. The solution is identical to the Earth-to-Moon 
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solution of the previous section. Assuming the velocity, v(t^), just prior to 
the midcourse maneuver is known, the velocity correction is given by 


Av(t o ) 




^ t o ) 


(B3-1) 


where 'v(t ) is determined from the asymptotic solution. 

I 

B3. 3 Interplanetary Midcourse 

The interplanetary midcourse problem is shown in Figure B8. The initial 
time, t , the initial position relative to the sun, r^t Q ), and the pericenter 
radius, p^, inclination, i^,, and time, t ,p f at the target planet are all pre- 
scribed. The initial velocity relative to the sun, v(t ), is unknown and must 
be determined. 


This problem is also solved using the midpoint-to-target body solution of 
Section B2. 1. The primary body is the sun and an ephemeris is required 
giving the planetary motions in cartesian coordinates with the sun at the 
origin. In this problem the subscript T refers to the target planet and the 
effects of all other planets enter only through the constants 'V_ r p, 6_,p, 
and q^. 


B 3 . 4 Inte rplanetary 

The inte rplanetary boundary value problem is shown in Figure B9. The peri- 
center radius, p^, inclination, i^, and time, t are prescribed at both the 
launch planet, k = L, and at the target planet, k = T. The hyperbolic excess 
velocities at both planets are unknown and must be determined. 


This problem is solved using the launch body-to-target body solution of Sec- 
tion B2. 2. The primary body is the sun and the same ephemeris as is used 
in the interplanetary midcourse solution is required. In this solution the 
perturbing effects of the planets from t = t^ to t = t Q are included in the con- 
stants V T , 6 T , and q T while the effects from t = t to t = t_ are included 

— J_/ i—j OX 

inY_ T , 6_ t , L t and q_ T . 
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Figure B9. Interplanetary Solution 


B4 SPECIAL MOON-TO-EARTH SOLUTIONS 

In this section two special moon-to -earth problems are considered. In both 

problems the boundary conditions are the initial and final times, t. and t , 

I e 

the initial position relative to the moon, R(t^), the entry radius relative to the 
earth, r^, the inclination relative to the earth, i , and the entry flight path 
angle, Y e . Also in each problem the trajectory prior to t = t^ is assumed to 
be an orbit about the moon with velocity V(tJ) just prior to t^. The first prob- 
lem involves finding the velocity V(t^) which results in a trajectory satisfying 
the earth entry conditions. This constitutes a single impulse problem where 

AVj = V(t+) - V(t“) (B4- 1 ) 

is the impulsive velocity. 

In the second problem the velocity after the impulse is assumed to be of the 
form 

V(tJ) = (1 + Ij) V(t“) (B4-2) 


which gives an impulsive velocity of 

AVj = Ij V(t“) (B4-3) 

The new velocity does not necessarily result in a trajectory satisfying the 
earth entry conditions and a second impulse AYL, is applied at t = t^ where 
(t^ - tj) is small compared to the total flight time, (t - t^). The second 
impulse must give a trajectory satisfying the entry conditions. This consti- 
tutes the two impulse problem. 

B4. 1 Modified Lambert Problem 

In the standard Lambert problem two position vectors are prescribed as well 
as the flight time from one position to the other. Solution of the problem 
gives the two-body trajectory connecting the two positions. Lambert's 
theorem is stated in functional form as 
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(B4-4) 


t 2 - t 1 = F (a, x 1 + x 2> c) 

cfi- 

whe re 

2 2 1/2 

c = (Xj + x 2 + 2 xjX 2 cos 9) (B4-5) 

■*' ... -• .! i. 

When the two position vectors are given 

cos 0 = (xj • x 2 )/(xjX 2 ) (B4-6) 


Solution of the Lambert problem requires solving (B4-4) for a. Once a is 
known a set of equations is solved to eventually give x(t). Several forms of 
(B4-4) have been proposed, all of which require numerical techniques to 
solve for a. 


For certain problems with prescribed entry conditions the final position x., is 
not known, only its magnitude is given. However the inclination and entry 
flight path angle are also prescribed and allow for a solution. The angle 0 
between the two positions is defined as 

0 = f 2 - fj (B4-7) 

where f^ and f 2 are the initial and final values of true anomaly given by 


fj = cos 


-1 


f 2 = cos 



. , TT < f^ £ 2 tt 


, TT £ f 2 2 2 tt 


(B4-8) 


(B4-9) 


The range of f^ and f 2 comes from restricting the flight time to be less than 
that for 0 = tt and having the entry before perigee. The eccentricity can be 
written as 
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j^a 2 + (x^ - 2a) cos 2 Y 2 J * ^ j ^ 


(B4-10) 


where Y^ is the flight path angle at t^. This modified Lambert problem 
requires a simultaneous solution of (B4-4), (B4-5) and (B4-7) - (B4-10) for 
a, c, 0, f^, f^ and e, i.e., six equations for six unknowns. 


.. 5 


In Section B2. 1 the approach of Battin was used to determine the ascending 
node of the inner solution. A similar approach is now used to determine the 
ascending node of the modified Lambert solution. Let the initial position x^ 
have components (^, q ^ ) and, using Figure B2 with replaced by x^, 


define 


. , /, t 2 , 2 . 1/2 

sin X = nj/Uj + nj) 


:os\ = ij/(£i + ni ) 1/2 


tan o = ^/(£ 2 + q 2 ) 1/2 


(B4- 1 1) 
(B4-12) 
(B4-13) 


Since the plane of motion must contain x^ the ascending node must satisfy 


Si = X. ± cr + ( 1 ± 1 )tt/2 


(B4-14) 


where 


cr = tan a /tan i 


(B4-1.5) 


and i is the prescribed inclination. The ± sign indicates two possible solutions 

5 

satisfying the prescribed inclination. The unit normal is now defined by 


n = (sin £2 sin i, -cos Si sin i, cos i) 


(B4-16) 


Using standard formulas for elliptic motion (c-f Battin) the initial position and 
velocity are 


xftj) = Xj 


(B4-17) 
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where 


V (tj) = e sin fj,/(nab) 

V (t.) = (1 + e cos f.)/(nab) 

<p l i 

and 


n = a 


-3/2 


b = a( 1 - e 


2 , 1/2 


The final position and velocity are 


x(t 2 ) = x 2 


X 2 

= (cos 0 Xj + sin 9 n x x^) 


i (t 2> = ST [ V r<y S2 + V ^ (t 2> 2- X 5 e] 


where 


(B4-19) 

(B4-20) 

(B4-21) 

(B4-22) 

(B4-23) 

(B4-24) 


V (t 2 ) = e sin f 2 /(nab) 

V ,(t-) = (1 + e cos f-,)/(nab) 

<p C. c 

The solution at any time t is given by (A6-2), (A6-5) and (A6-?6), i.e,. 


(B4-25) 

(B4-26) 


x(t) = f(t) xjtj) + g(t) kitp 


(B4-27) 



whe re 


f(t) = 1 - a(l - cos AE)/xj (B4-28) 

g(t) = (x^ sinAE)/(na) + e sin (1 - cosAE)/n (B4-29) 

and 


AE = E - E, 


(B4-30) 


The eccentric anomaly E is defined from 


n(t-t) = E-esinE 

P 


(B4-31) 


or by 


sin E = (1 - e^)*^ sin f / ( 1 + e cos f) 


cos E = (cos f + e)/(l + e cos f) 


(B4-32) 

(B4-33) 


Since e and f^ are known (B4-32) and (B4-33) can be used to determine E^. 
Then from (B4-31) the time of pericenter passage, t^, is 


tj - (E^ - e sin E^)/n 


(B4-34) 


The equations derived in this section completely define a two-body ellipse 
satisfying the prescribed values of the initial position, x^, the inclination, 
i, the final radius, the final flight path angle, V and the time t^ - t^, 

where both t^ and t^ are given. 

B4. 2 Single Impulse Solution 

B4. 2. 1 Outer Solution 

The boundary conditions which are prescribed for the outer solution are the 
earth entry time, t e> radius, r e> inclination, i^, and flight path angle, 
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Thus in the modified Lambert problem the subscript 2 is replaced by e, i.e., 


x_ = x 


i-, = i 


V, = V 


(B4-35) 

(B4-36) 

(B4-37) 

(B4-38) 


An additional boundary condition for the outer solution is the initial time t.. 
The initial position for the modified Lambert problem is taken as the position 
of the moon at t = t . , i.e., 


-1 = £m ( V (B4-39) 

This forces the Lambert solution to pass through the center of the moon as in 
the earth-to-moon solutions. ' 


Solution of the modified Lambert problem gives the zeroth order outer solution 


from (B4-27), i. e. , 

r^(t) = x(t) (B4-40) 

The position and velocity at t^ are given by (B4-17) and (B4-18) 

r^tj) = x(tj) (B4-41) 

^(tj) = x(tj) (B4-42) 

and at t by (B4-23) and (B4-24) 

= x(t 2 ) (B4-43) 

^(tg) - x(t 2 ). (B4-44) 
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The boundary conditions on the higher order outer solutions come from (B2-29) 


and (B2-30) 

£(t e ) = £ 0 (t e ) + ^(tg) + i± Z £ 2 (t e ) 

= Zo^e') + M-£i(t e ) + p 2 £ 2 (t e ) 


(B4-45) 

(B4-46) 


The effect of the higher order boundary conditions will be considered shortly. 
No further information about the outer solution is required at this point. 


B4. 2.2 Inner Solution 

The only boundary condition prescribed for the inner solution is the initial posi- 
tion R, , at t = t,. However R, , is actually a function of the inner variable 
— M 1 — M 1 


S M (t " ‘pMV^ 


(B4-47) 


where t ^ -is the time of pericenter passage of the inner solution and p= 
the dimensionless mass of the moon. From (A13-2) 


pM 


= t A/f + I1T . 

M ' M 


(B4-48) 


where t. , and r, , must be determined. (In the previous boundary value solu- 
M M 

tions it was stated that t could be set equal to zero without loss of generality. 

The effect of a non-zero t will be demonstrated in the next section. ) When 

t = t, then S., = S.,, and the initial position is R(S W1 ). 

1 . M Ml r - Mr 

From (A12-54) the hyperbolic excess velocity is 


V = ( U u> V A/r> 

— ooM ' M M M 


(B4-49) 


The excess velocity and the initial position are sufficient to determine the 
inner solution as follows: 


140 



The position and velocity at any time S are given by 


R. .(S..J = R»,(Sw! cos f. .(S.. .) e + R, , (S. . ) sin f (S ) e, 

— M x M M M M M —a M M M M — b 




M 


sin f, ,(S, .) e + e w + cos £. ,(S W ) 
M' M — a 1_ M M' M 


K] 


where 

r ’ • 

II 

| — M 


[*M < 


1 / 211/2 


(B4-51) 


(B4-52) 


(B4-53) 
(B4-54) 

and f^. is the true anomaly. The unit vectors e^ and e^ lie in the orbital plane 
with e directed toward pericynthion as shown in Figure BIO. 

ct 



Figure BIO. Orbital Plane Coordinates for Inner Solution 
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Now let 


-Ml = -M (S M1 
-Ml = -M (S M1 
*M1 = f M (S Ml* 


) 

) 


Then (B4-51) and (B4-52) may be solved for e and e, giving 

“3. — D 


e M + cos f Ml , „ 

\ -A I -Ml 

M 


R Mlf in ( M I 


t-M 1 


% 


'sin f... \ /R ul cos f ul 

ML n / -Ml. MI . 

-2 ) — M 1 \ 7., 1M1 


M 


M 


Substituting (B4-58) and (B4-59) into (B4-51) gives 

5m = [‘ 


2 - Hu + 'OS (f M - f M1 )] R M1 /< 2 


M M 


M 


+ R M R M1 Sin (f M ' f Ml' -Ml /f M 


Now define e by 

— 00 


lim / — M 

= r m“"\ r m, 


and <p by 


^ " R - co (f M " f Ml' 
M 


(B4-55) 

(B4-56) 

(B4-57) 


(B4-58) 


(B4-59) 


(B4-60) 


(B4-6 1 ) 


i 


(B4-62) 
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Then, from (B4-60) and (B4-61) 


1 ^ „ . Ml . . 

7T (1 - cos 0) ^11 + "T~ SIn * Y.M1 
'm M 


(B4-63) 


Solving for V. ^ gives 


-Ml 


R 


g M -coM (1 - cos 0) -Ml 

Mi sin 0 V “m 7 M sin ^ R M1 


(B4-64) 


where e has been replaced by V , r using 
— ® r 1 — coM b 


e ' = V /V .. 

— oo — coM CO M 


(B4-65) 


i. e. , the radius vector becomes parallel to the velocity vector at infinity 
(cf. (A12-58)). 

Now define Y^ and the unit vectors e^ and e_^ as shown in Figure Bll. 

The velocity V. , , can be written 
—Ml 


V,,. = V... cos Y. e + V,,, sin Y.'e. 

—Ml Ml 1 —a Ml 1 — p 


(B4-66) 


Also 


e = cos be + sin 0 e r 

— CO r — Q r —0 


(B4-67) 


Eliminating e_^ gives 


V A ,, = V,,, (cos Y, - sin Y ,/tan 0) e +V W1 sin Y,/sin 0 e (B4-68) 

— M 1 Ml 1 1 —a Ml 1 r — ® 


Using (B4-65) and 


e = R. , /R 
— ot — M 1 Ml 


(B4-69) 
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Figure Bl 1 . Orbital Plane Coordinates at t = , 

give s' 

/ sin R-jl r i sin V V 

V = V ( cos Y - ) — + V i 

~ M Ml V 1 tan ^/ R M1 M1 - sin V V « M 

Equating coefficients of l n (B4-64) and (B4-70) gives 


(B4-70) 



R 


Ml V M1 


sm 


(B4-71) 


t • 

which is simply the definition of the .scalar angular momentum which agrees 
with (B4-54). Equating the coefficients of Rj^ gives 


sin Y 


V.,, cos Y. - — 

M 1 \ 1 tan (p 


1 


( 1 - COS 0) 

I M tan. <p 


(B4-75) 
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Eliminating using (B4-71) and solving for gives 


-L 

M 


Ml 


R 


M COS 0 (1 - COS 0) 

sin <P 


R., . L Ml si n 0 M 
Ml 


(B4-76) 


From the conservation of energy 


Ml 1 


V‘ 


.M 


R 


Ml 


(B4-77 ) 


Substituting (B4-76) into (B4-77) and solving for gives 


r mi v <x,m sin * 


M 


1 + /I + 


R 


Ml V coM (1 +COS \ 


(B4-78) 


where the plus sign is used in front of the radical in order to satisfy (B4-54). 


Using (B4-78) in (B4-64) eventually gives 


-Ml 


J_ 

2 


1 + /I + 


r mi v »m (1 + cos V 


V A/r 
— coM 


'M 


2 R 


Ml 


1 - /l + 


?M1 V .M (1 + cos *• 


R» , , 
—Ml 


(B4-79) 


This expression gives the initial velocity V... in terms of V . , and R, ,, since 

c ’ — M 1 — ooM — M 1 

i 

cos * - <«M1- ' < R M1 V „M> ■ (B4 - 80) 
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Finally the inclination is defined by 


C0S *M = 2.3 ’ Sm 


(B4-81) 


where is the unit normal to the orbit plane and is defined by 


-M " ( -Ml X -Ml* / 


(B4-82) 


Therefore, combining (B4-81) and (B4-82) gives 


cos = <2 3 • —Ml x ^M1 ,/( M 


<B4-83) 


Substituting (B2-11) and (B4-54) into (B2-10) gives 





, 1/2 



(B4-84) 


It was shown in Section B2 that prescribed values of the hyperbolic excess 
velocity, the pericenter radius and the inclination are sufficient to determine 
the inner hyperbola and its two constants, A.^ o and C^ o . Equations (B4-83) 
and (B4-84) give the radius and inclination as functions of the excess velocity 
and the initial position vector. It has therefore been shown that the excess 
velocity and initial position do indeed determine the inner hyperbola as was ' 
stated at the beginning of this section. 5 


A final parameter which must be determined before proceeding to the boundary 
value solution is t^.. In the previous boundary value solutions it was stated 
that was arbitrary. Such is not the case now, as can easily be shown. The 
hyperbolic eccentric anomaly, F, at time S^. = is given by the standard 

expressions (cf. Carlson^). 


cosh F 1 = (a M + R Ml )/(a M e M> 


sinh F i = i^r a M 1/2) 


(B4-85) 

(B4-86) 
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Then from Kepler's equation, (A9-2), 


S M1 


(e M sinh F, - Fj) / n M 


and, from (B4-47) and (B4-48) 

(S S.. , , = (t. - t w - u w t w ) / p., 

' Ml 1 M r M M 


(B4-87) 


(B4-88) 


The time t. , is the time at which the zeroth order outer solution passes 
M 

through the center of the moon. But from (B4-39) 


t 


M 


t 


1 


and therefore, from (B4-88) 


(B4-89) 


T M ~ " S M1 


(B4-90) 


(Note: It would be possible to choose t^ different from tj and. still obtain a 
solution. However, the value of is required before can be calculated 
because the inner solution requires knowledge of the zeroth order outer solu- 
tion to get the zeroth order approximation for This means that tj, 

and Sj are always fixed in (B4-88) and is then given by 


M 


= - S 


Ml 


< fc l - 


t M ) /|a 


(B4-91) 


;, Thus is never arbitrary. (The choice of t^ = t^ gives the simplest result. ) 
B4. 2.3 Boundary Value Solution 

The single impulse boundary value problem is shown in Figure B12. It is 
^'similar to the earth -to -moon problem shown in Figure B6. The terminal 
conditions at the earth are given by (B2-27) and (B2-28) with t Q = t . From 
/•Figure B 12 

r^tg) = r(t e ) , (B4-92) 
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therefore 


^(tg) = £ 2 (t e ) = 0 


(B4-93) 


and, as in (B2-3 1), let 


8 v(t e ) = Xi^e^ + f 1 


(B4-94) 


Now let k = M with 


H- = H 


M 


(B4-95) 


Q M = +1 


(B4-96) 


and substitute (B2-24), (B2-25), (B4-93) and (B4-94) into (Bl-8). Since 


t^ = t^ (Bl-8) reduces to 


B^.y 8 v(t e ). = L m 


1 




T , - + log -pr 

M n M \ 2n M 


— coM + — M + (B4 " 97) 


-1 


D(t 1 ,t e ) 8v(t e ) = n (V a)M -V M ) + 6 M + ,_ M 


(B4-98) 


where 


— M 


= ^o (t i> "£m (t i } 


(B4-99) 


should not be confused with defined by (B4-52). Solving for 6v(t g ) 

and V ,, gives 

— coM G . • 


6v(t e ) = B(t r t e ) 


-1 


M 


.*1", / ^M e M 
T M h M og \ 2n M ! 


V Tv r 
— coM 


+ Xm + * £m 


(B4-100) 


149 



(B4-101) 


-co M ” -M + ^ D(t l’ te) 6 - (t e ) " ^Am " R 3-M 

These two expressions are not explicit since V ^ and 6 v(t g ) appear on the 
right hand sides. The solution is obtained from the following steps: 

1. Given the initial time t^ solve (B4-39) and (B4-41) for r^tj) using 
the ephemeris for 

2. Using the prescribed terminal values of t , r , Y and i solve the 

c c c c 

modified Lambert problem given in Section B4. 1. This gives v^t^), 

r (t ) and v (t ). 

— o e — o e 

3. Solve (B4-99) for where obtained from the ephemeris. 

■r 

4. Let V . , = V. This gives the zeroth order excess velocity. 

— coivl — 1\A 

5. Using the prescribed initial position R ^ evaluate (B4-78) - (B4-84). 

6. Evaluate (B2-1) - (B2-26) and (B4-85) - (B4-90). 

7. Evaluate (B4-100) using the results of step 6 and with = 0. This 

gives the first order 6v(t g ). 

8. Evaluate (B4-101) using the results of step 7 and with = 0. This 

gives the first order excess velocity. 

9. Repeat steps 5 and 6 using the first order excess velocity. 

10. Evaluate (B4-100) using the results of step 9. This gives the second 
order 6v(t e ). 

11. Evaluate (B4-101) using the results of step 10. This gives the 
second order excess velocity. 

12. Repeat step 5. This gives the second order initial velocity V j. 
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(B4-102) 


13. The single impulse is given by 


AV W1 = V... - V„(S A/1 ) 
—Ml —Ml — M' Ml 


where (S^. ^ ) is the velocity just prior to the impulse. 


14. The entry position and velocities are given by (B4-92) and 


v(tj = v (t e ) + u6v (t ) 


(B4-103) 


Steps 12, 13 and 14 represent the solution of the single impulse problem using 
the standard technique. A nonlinear solution is also possible using 


Eo (t l> = lo'V + ^ 


L u 
— M 


^M e M' 


+ KT7-K* 


M n 


M 


M / J 


V 

— coM 


+ — M + M 


(B4-104) 


instead of (B4-39) - (B4-41) as the endpoint for a new modified Lambert 
problem. Solution of the new 
Then (B4-101) is replaced by 


problem. Solution of the new Lambert proglem gives v^(t^) and v^(t e ). 


-coM -M “ M ' ^ — M 


(B4-105) 


where 


-M 


^(‘l’ 'Em (t i! 


(B4-106) 


and (B4-103) is replaced by 




(B4-107) 
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The non-linear solution must also be evaluated in the sequence of the standard 
solution since the right hand side of (B4-104) is a function of . In practi- 

cal applications the entry will always be near perigee and the non-linear solu- 
tion is probably preferable to the standard solution since it bypasses the 
errors introduced by the linear matrix B(tj, t^). An even greater advantage 
of the non-linear solution over the standard solution is the fact that the latter 
does not satisfy the entry conditions exactly while the former does. 

In the linear solution, the modified Lambert solution (the zeroth order solu- 
tion) satisfies the entry conditions exactly but because of (B4-94) the first and 
second order entry velocities are slightly different. Thus the entry conditions 
remain satisfied only to order unity. In the non-linear solution the total entry 
velocity, 

v(tg) = v^tg) + p Vj(t e ) + p 2 v 2 (t g ) (B4-108) 

is replaced by (B4-107) which comes from the modified Lambert solution. 

Since the Lambert solution satisfies the entry conditions exactly so does the 
total entry velocity and the entry conditions are satisfied to second order. 

The non-linear solution is shown in Figure B13. 

B4. 3 Two Impulse Solution 

The two impulse problem might be thought of as an extension of the single 
impulse problem. It is, however, a somewhat different problem since the 
impulse which gives a trajectory satisfying the earth entry conditions occurs 
not in the inner region, as in the single impulse case, but rather in a region 
which is more closely associated with the overlap domain. Since both the 
inner and outer solutions are valid in the overlap domain it should be possible 
to represent a trajectory in this region by either solution. This leads to two 
alternate approaches to solving the two-impulse problem. These approaches 
are developed in the following sections. 

B4. 3. 1 Inner Solution 

The inner velocity prior to the first impulse is Then according to 

(B4-3) the first impulse is 


152 




153 



h — M < S I> 


(B4-109) 


AVi 


The value of 1^ must be chosen to give a hyperbolic trajectory after the 
impulse, thus 


< s ? - 1 


(B4-110) 


where is the initial inner radius. Then the velocity following the 

impulse is 




(1 


+ I 1> 


— M 


(Sj) 


(B4- 111) 


The velocity (S^), which is assumed known, uniquely defines the zeroth 
order inner hyperbola. The elements defining this trajectory can now be 
obtained. They are 


i = < S 1> X < S j> 

= angular momentum vector 



= angular momentum 

n = i/T 

= unit normal vector 
h = V M (sj ) 2 /2 - 1/R m (Sj) 
= energy 
a = 1 /(2k) 

= semi-major axis 


(B4-112) 


(B4-113) 


(B4-114) 


(B4-115) 


(B4-116) 
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= mean motion 
i = cos” * (N • £ 3 ) 

= inclination 

In addition, there are the argument of the ascending node, 
sin f 2 = (N ' £j) J sin i 
cos = - (N * £^) / s ^ n i 

the initial eccentric anomaly, 

sink F [ = [r m (S,)- V M (st)]/(I 1 / 2 e) 

cosh F x [ a + R m (Sj) J/ a e 
the initial true anomaly, 

sin Tj = ^a (e " 2 - l ) 1 ^ 2 sinh R M (S^ 

cos fj = a(e - cosh Fj) j R^ (S 
and the argument of pericynthion 


(B4-117) 

(B4-118) 

(B4-119) 

(B4-120) 

(B4-121) 

(B4-122) 

(B4-123) 

(B4-124) 

(B4-125) 

(B4- 126) 
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where 


sin co 


|W*NxR m (S,)] 
R(S ^ ) sin i 


- g - (S 1> 

w ~ R(S ) sin i 


(B4- 127a) 


(B4- 1 27b) 


All of these expressions are derived from standard two-body relationships. 
Finally, from Kepler's equation for a hyperbola 


= (e sinh - F ) / n 


while the inner time is defined by (A7-27) and (A 13- 2) as 



(B4- 128) 


(B4-129) 


Equating (B4-128) and (B4-129) gives 

t M = (t l - V 7 V- + (F 1 - * sinh F l ) (B4-131) 

In this particular problem t^ is fixed and t M is arbitrary. Putting fc M = tj 
gives 

t., = (F. - e sinh F.) / n (B4-131) 

Ml 1 

Using the elements defined here it is possible to determine the behavior of 
the inner solution when S is large using the results of Section A 12. The 
values of A', B' and C' are found from (A12-27) - (A12-29), the values of 
U' and V' are found from (A12-38) and (A 12-39) with Q = 1, and the values 
of V .... and L,. are found from (A 12-48) - (A 12-55). 

— co M — M 
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The values of V ,, and can be used in the fundamental solution to deter- 
— coM — M 

mine the position and velocity at some later time, i. e. , the time of the second 
impulse. This approach would be followed if the second impulse were to be 
added strictly in the outer domain. 

If the second impulse occurs in a time period approximately 1 1 to 18 hours 
after the first impulse then it occurs in the overlap domain where both the 
outer and inner solutions are valid. Therefore, rather than using the funda- 
mental solution, which includes the matching of the inner and outer solutions, 
the position and velocity prior to the second impulse can be determined from 
an inner, perturbed hyperbola alone. Carlson has derived some' simple 
formulas to achieve this. 

The position and velocity at any time S, where S > 1, can be found from 


5*1 (S) = 5 Mo (S) + H 2 R Mz (S) + ^ R Ms (S) (B4-132) 

5m < S > = 5 Mo <S) + / Y^S) + ^ V M3 (S) (B4-133) 

where and —Mo are t * ie t w °-body position and velocity, 

5m 2 < s > = g m[5mo< s)s2 /2-5mo ( s » s3 / 3 ] (B4 - 134) 

5m 2 < s > = g m [5 Mo (S) s - Ymo (s > s2 / 2 ] (B4 - 135) 


and 



a 3 s 


(B4- 136) 



4*3 S 3 


(B4- 137) 
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The two-body position can be obtained from any standard expression, (A9-1) 
is an example. The two-body velocity can then be obtained by differentiation 
with respect to S or by standard velocity expressions. The equations for 

an d Y.M2 are s ^ m P^^ cat ^ oris of Carlson's expressions found by substitut- 
ing G^, which is defined by 

= G( £M (t l )) (134- 138 ) 

for his time averaged gravity gradient matrix. The expressions for amcl 

—M3’ ^ ounc ^ ^om (A12-102) and (A12-103), make a first order correction for 
the substitution of the constant matrix for the time averaged value which 
he used. Using G^. simplifies his expressions considerably with very little 
loss of accuracy. 

The expressions for Rjy^ an< ^ —M2 were derived using Taylor series expan- 
sions^ while the expansion for R^^ Section A12 is an asymptotic expansion. 
The two forms are equivalent if the product nS » 1, however the algebra neces 
sary to show this is long and tedious. For moderate values of nS, (A 12-81) is 
not a good approximation and (B4-134) is preferred. , 

* 

B4. 3. 2 Overlap Solutions 

The second impulse is to be applied at some prescribed time t = t^. The 
impulse must result in a trajectory satisfying terminal constraints at a 
later time. In order to calculate the impulse it is necessary to determine 
the position and velocity at t^, i. e., at the instant just prior to adding the 
impulse. There are two possible ways of achieving this. 

The first method involves the use of the fundamental solution. This requires 
knowledge of a zeroth order outer solution. As in the single impulse solution, 
(B4-39) - (B4-41), let 

W = £m <*i> (B4 - 13,) 
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This gives the initial position. Now solve (B4-99) giving 


v 

— o 


{t l ) “ (t l* + -M 


(B4- 140) 


This is the initial velocity. In other solutions V^. has been used as the first 
approximation to Since was determined in the previous section 


let 



Y 


coM 


(B4- 141) 


Then at any time t the zeroth order solution is given by (A6-2) and its 
derivative 


r Q (t) = f Q (t) r Q (t 1 ) + g Q (t) V q (t x ) (B4-142) 

v (t) =..f (t) r (t ) + g (t) v (t ) (B4-143) 

— o o — o 1 o — O 1 

These expressions can be used to find the zeroth order position and velocity 
at t = t^- 


Since the moon corresponds to the launch body the perturbations at t^ are 
found from (B2-56) by letting t = t^, L = M and adding the second order 
terms from (Bl-8). The result is 


— 1 ( t 2 ^ + ^ -2 ^ 2 ) 


1 (t 2 } + ^2 (t 2 } > 


where the constants V.„ 

— M 



— M’ — M atK ^ — M are eva ^ ua *- ec ^ along r^(t). 
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The position and velocity at t = t^ are found by combining (B4- 142) - (B4-144) 
giving 

r(t“) = L q ^~ 2 ) + + ^ 2 £ 2 (t 2 } (B4- 145 ) 

v(Q = V o (t-) + ,v l( t-) + , 2 X2 (t-) (B4-146) 

By considering all of the error terms in Sections A12 and A 13 it can be shown 
that in the region of interest the dominant error in (B4-145) is O^j, / (t^ - tj) 2 j. 
The following table for position error can be constructed: 


t^ - t . , hours 

5 

11 

18 

20 

0(t 2 - tj) 

2/3 

1/2 

2/5 ‘ 

3/8 

8 - 

M- 

H- 

M- 

0(,x 3 /(t 2 - tj) 2 ) 

1. 67 

2. 00 

2. 20 

2. 25 

M- 

M- 

M- 

M- 


Thus as t 2 - tj increases the theoretical error decreases (i.e., the exponent 
of p increases). Although the first term ignored in (B4-145) is order p^ the 
error shown in the table is somewhat larger. 

The second method does not use the fundamental solution but rather the per- 
turbed hyperbola discussed in the previous section. Solving (A7-26) for r_ 
gives 

— = -&M + M (B4-147) 


and differentiating with respect to t gives 


d^r 

dt 


dR 




M 


dS 


(B4- 148) 
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(B4- 149) 


Letting chr/ dt = _v(t) and dR^/dS .= V^(S) gives 

v(t") = £ M (t 2 ) + -M (S 2 ) (B4-150) 

where and are the position and velocity of the moon at t = t^ and 
R^-CS^) and are obtained from (B4-132) - (B4-138). The value of 

is found from 

s 2 = s i + (t 2 - t L ) / p (B4-151) 

Based on Carlson's results^ the dominant error in (B4-149) is 0^p(t 2 - t^) 2 ^ 
which results in the following table for position error: 

t^ - t^, hours 
0 (p(t 2 - t^ 2 ) 

This solution has an increasing theoretical error as t^ - t^ increases (i. e. , 
the exponent of p decreases), which is just the opposite of the fundamental 
solution. At 5 hours the perturbed hyperbola is more accurate, at 1 1 hours 
the accuracy is the same for both solutions, and at 18 and 20 hours the funda- 
mental solution is more accurate. It would therefore appear as if one solution 
or the other should be used depending on the value of t^ - t^. The following 
points, however, should be, considered: 


11 


18 


2.33 2.00 1.80 

L P- ^ P 


20 

1. 75 


1. The perturbed hyperbola solution was evaluated for a hyperbola of 

2 81 

moderate energy. At 5 hours the actual error was p ’ rather 

2 33 2 05 

than the predicted p * . At 20 hours the actual error was p * 

1 75 

rather than the predicted p ' . Thus the actual errors are less 

than predicted. 
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2. Numerical analysis has shown that errors associated with the second 
order fundamental solution are usually larger than predicted. 

3. The perturbed hyperbola solution is much easier and less time con- 
suming to evaluate since it does not involve any definite integrals and 
uses basically two-body functions. 

Although either (B4-145) and (B4-146) or (B4-149) and (B4-150) give the posi- 
tion and velocity at t^ the latter is preferable. 

B4. 3. 3 Outer Solution 

The outer solution will be that part of the trajectory which follows the second 
impulse. The boundary conditions are the position r^t^) and the entry condi- 
tions of time, t , radius, r , inclination, i , and flight path angle, V . In 
6 6 C 6 

the modified Lambert problem the subscript 1 is replaced by 2 and the sub- 
script 2 replaced by e so that 



(B4-152) 
(B4- 153) 
(B4- 154) 
(B4- 155) 
(B4- 156) 
(B4-158) 


Solution of the modified Lambert problem gives the zeroth order outer solution 

r (t) for t > t_. 

— o 2 

The first and second order outer solutions are obtained directly from (A6-10) 
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giving 


'r^t) + fir £ ( t ) 

^(t) + (J.v 2 (t) 




£ L ( y + nz 2 (y 
,v,(t e i + ^ 2 <t e y 



where 



(t, t £ ) 


-11 


(t, t £ ) 


r 20 (t « fc e > 


-21 


(t, t £ ) 



B(t, t) F x (t) dr 


D(t, r) F L (t) dr 


B (t, T ) F 2 (t) dT 


D(t, T) F Z (T) dT 


B4. 3. 4 Boundary Value Solution 

The two-impulse boundary value solution is shown in Figure B14. 
nal conditions at the earth are given by (B2-27) and (B2-28) with t 
Figure B14 

r (t ) = r(t ) 

— o e — e 


(B4-159) 


(B4-160) 


(B4-161) 


(B4-162) 


(B4-163) 


The termi- 

= t . From 
e 

(B4-164) 
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0 ) 

co 



(- 1 1 
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so that 


Ii(t e ) = *2 (t e > = 0 


(B4-165) 


Since the velocity perturbations do not vanish let 
6v(t e ) = Vj(t e ) + pv 2 (t e ) 


(B4-166) 


+ . 


The outer position and velocity at t = t is then 



(B4-167) 


However, from Figure B14 



v(t 2 ) + AV 2 


(B4-168) 


and 


i (t 2> = * £ (t 2> 


(B4- 170) 


where t ^ ^ * s ^ our id from the inner solution. Then (B4-167) reduces to 


6v(t e ) - -B(tJ, g- 


[-10^2’ + ^-20^2* fc ej] 


(B4-171) 
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(B4-172) 


^ — 2 = Io (t 2> - 2(‘ 2 ) + I* [mtj, t e ) 6v(t e ) 

+ A i ( l 2» l e *] + ^ ^2 1 * l 2' *> * 

The solution is obtained from the following steps: 

1. Given the initial position R^.(Sj), velocity V^(S^) and impulse 1^ 
evaluate (B4-111) through (B4-131). 

2. Using the results of step 1 and the appropriate equations in Section 

A 12 find the position and velocity at t^ from either (B4-142) - (B4-146) 
or (B4- 149) - (B4-151). 

3. Using (B4-170)for the initial position and (B4-154) - (B4-157)for the 

terminal conditions solve the modified Lambert problem of Section 

B4. 1. This gives v (t 0 ), r (t ) and v (t ). 

' e — o 2 — o e — o. e 

4. Evaluate (B4-171). This gives the correction to the entry velocity. 

5. Evaluate (B4-172). This gives the second order AV,,, ^he secon d 
impulse. 

6. The entry position is given by (B4- 164) and the entry velocity by 

v(t e ) = V Q (t e ) + p6v(t e ) (B4-173) 

Steps 5 and 6 represent the standard solution for the two impulse problem. As 

in the single impulse solution the entry conditions are not satisfied exactly 

because of the non-zero 6v(t ). A non-linear solution which satisfies the 

— e 

entry conditions exactly to any order is obtained by first rewriting (B4-167) 
and (B4- 1 68 ) 
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d‘ 2 > 


X (t 2>‘ 


+ 


~ M- 


A ^2, 


— 1 0-^2’ fc e } 


,-ll (t 2’ t e ) 




—20^2* t e ) 


k -21 (t 2’ t e ) 


r (tt) 

— O L 




+ n®(t 2 , t ) 


(B4-174) 


kSv(t e )> 


The right hand side of (B4-174) is the sum of a two-body solution plus the 
propagation of initial variations at t Thus the right side is a pure two body 
function and can be replaced by 



0 


AV 
-2 


- h • 


—10^2* fc e* 


,^ll (t 2’ t e ) 


- h 


—20^2’ t e ) 


—2 1 ^2’ t e ) 


^ (t 2 } 


—o^2^' 


(B4-175) 


This defines a new modified Lambert problem for which the new initial posi- 
tion replacing (B4-170) is 

x 'l = r o (t ? ) 

— l — o 2 


— ( fc 2? " fJ '-10 (t 2’ fc e* " ^ —20^2’ t e ) 


(B4-176) 


Solving the new modified Lambert problem between x' and the entry conditions 

i + ? - 1 

gives v (t ^ )» (t ) and v' (t ). From (B4-175) the velocity impulse is 


A — 2 ~ -o^2^ “ -^2^ + ^11^2’ t e ) + ^ ^21^2’ t e ) 


(B4-177) 
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while the entry position and velocity are 


r (t ) = r ' (t ) (B4-178) 

— e — o e ’ 

v(t ) = v' (t ) B4-179) 

— e — o e 

Since (B4-176) and (B4-177) come from the modified Lambert solution they 
will satisfy the entry conditions exactly to any order! The non-linear solution 
is shown in Figure B15. 



Figure B15. Non-Linear Version of Two-Impulse Moon-to-Earth Solution 
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Section C 

EVALUATION OF PERTURBATION TERMS 
Cl TYPES OF PERTURBATION TERMS 

The n order outer solution is given by (A6-10). It consists of two types of 

til 

terms: (1) two-body propagation of the n order initial conditions, and (2) 

til 

integrated effects of the n order perturbations due to the N-2 perturbing 
bodies. Both of these effects enter into the boundary value solutions through 
the constants &_k> and Rk which in turn are functions of partial 

derivative matrices and definite integrals (plus algebraic and trigonometric 
terms). Formulas for evaluating the partial derivative matrices and the 
definite integrals are given in the following sections. 

C2 PARTIAL DERIVATIVE MATRICES 

th 

The partial derivative matrices are evaluated on the zero order outer 

solution r (t). - The formulas given here were derived by Goodyear 1 * and the 
~° 3 

notation is similar to that of Carlson. Two forms are presented; the first 
can be used for direct calculation of the matrices and the second is useful in 
evaluating the definite integrals. 

C2. 1 Goodyear Formulas 

The state transition matrix $(t , t) is a function of the two times, t and t. 

o o 

The eccentric anomaly difference . 

u = E q - E (C2-1) 

will be chosen as the independent variable. The Goodyear formulas then 
require the following functions of u: 

f( t o , t) = 1 - (1-cos u) (C2-2) 

o v o 7 
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2 . 

na sm u 

r (t ) r JV 
o o o 


g (t , t) = — r (t ) fa - (1-cos u) - sin u1 
6 o na o o L J 


e(t t) - 1 - a(1 ~ cos u) 

g ^o’ 1 r (t) 

o 


P(t Q , t) = 3(sin u-u) + e sin E q (1-cos u) 


+ (1-cos u) e cos E sin u 

o 


A RR (t o' t( 


r o 2( t) 


2 

(1-cos u) sin u 
r o (t) r o (t o ) 


A (t t) - (1 ~ cos u)sin u 
RV' o' ~ n r (t) r (t ) 
o o o 


a P(t ,t) . . . 

- n _ o ( 1 - cos u)sm u 

» t) — 


VR o’ 


n r 3 (t) n r 2 (t) 

o , o 


A ft t) - - (1 ' COS U) 
A VV o’ t ~ 2 ~~ 

n a r (t) 


B RR (t o’ fc) “ A RV (t o’ l) 


VV" - 

n a (t ) 

O O 


B VR* t o’ ^ A VV^o’ ^ 


(C2-3) 

(C2-4) 
(C2- 5) 

(C2-6) 

(C2-7) 

(C2-8) 

(C2-9) 

(C2-10) 

(C2-11) 

(C2-12) 

(C2-13) 
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1 


(C2-14) 


Byv^o’ *■) = — 3 ^ — 2 C^(sin u-u ) + { 1 -cos u) sin u] 
v J n a 


C RR (t o’ fc) 


na 


r (t) r (t ) 
o o o 


cos u _ , 1 

T 


r (t) r (t ) 2 ... 

o o o r (t) 
o 


1 


r 2( t ) 
o o 


sin u - 


na P(t Q , t) 

r 3 (t) r 3 (t ) 
o o o 


(C2-15) 


C Rv“o’ 11 


sin^ u ( 1 - cos u) 


2 ,. , r (t) r (t ) 
r (t ) L o o' o J 


(C2-16) 



0 0 


rt- 

0 . 

r"f~ 

II 

" A RR (t o’ fc) 

(C2-17) 

‘‘o’ 1 * = 

" A RV (t o’ fc) 

(C2-18) 

“o’ 1 * = 

C Rv' t o’ *> 

(C2-19) 

“o’ 1 * = 

( 1 -cos u)sin u a ^^ 0 ’ ^ 

(C2-20) 

nr o 2 “o* nr o 3 “o* 

(t 0 ,t* = 

“ A RV (t o' t| 

(C2-21) 

“o’ 1 ' =■ 

“ B RV (t o ,t) 

(C2-22) 


th 

where r (t ) and r (t) are the radii of the zero c order outer solution at t 
o o o o 

and t; a, e and n are the semi-major axis, eccentricity and mean motion of 
^(t); and 


a e . 

tr = — 77 — r sin E 
r (t o 

o' o 


*C2-23) 
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After some modification Goodyear's results can be expressed in the form of 
the transition matrix, i. e. , 


'A(t , t) B(t t) 
o o 


C(t ,t) D(t t) 

\ o o 


A*(t , t) B*(t , t)\ /E*(t , t) F* 
o o \ / o 


C*(t , t) D*(t , t) 
> o o 


F* E*(t , t) 
o / 


(C2-24) 


where 


A *( t o» t) - [§ (t 0 ’ fc) A RR (t 0 > t} A RV (t o’ fc) A VR (t o’ A W (t o’ t} ] (C2-25) 


B*(t o , t) 


■c*(t 6 , t)“ 


D*(t o , t) 


[" g (t 0 > B RR (t o’ ^ B RV (t o' ^ B VR (t o’ ^ B VV (t o! fc) ] 
[~K fc) C RR (t o’ fc) C RV (t o’ C VR (t o’ fc) C VV (t o’ fc) ] 
[ f(t o’ D RR (t o’ fc) D RV (t o’ D VR (t o’ fc) D .VV (t o’ fc) ] 


E*(t o , t) 




r (t ) v (t)' 
—o' o — o 


v (t ) r (t)' 
—o' o —o' 


v (t ) v (t)' 
—o o —o 


F* = 
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The elements of A*, B*, C*, and D* are all scalars while the elements of E* 

T 

are all 3x3 matrices, I being the identity matrix and xy being the outer 
product of x and (the superscript T indicating the transpose of the vector). 
Thus, each partial derivative matrix is the sum of five terms, each of which 
is the product of a scalar times a matrix. 


These formulas for the partial derivative matrices may be used either for a 

variable t, as in the definite integrals or for a fixed value of t. The unique 

characteristic of this formulation is that each partial derivative matrix can 

be written as a function of scalars and vectors (plus the identity matrix) which 

til 

are easily obtained from the zero order two-body solution. 

C2. 2 Modified Goodyear Formulas 

Although the independent variable u was introduced in the previous section 

the formulas still contained r (t), , v (t), and r (t). In the definite integrals 

— o o o 

it is advantageous to remove this dependence on t by introducing the modified 
derivative matrices^ 

r(t) 

B r (t o’ fc) = (C2-31) 


r (t) 

= — D(t o’ fc) 


(C2-32) 


and replacing r Q (t) and V Q (t) in the previous formulas by 


£o (,) = f < t o’ t) Io (t o ) + g( V t) io (, o ) < C2 - 33 > 

2o<‘> = f <‘o’ t) ^ (t o , + g(t o’ t) W < c2 - 34 > 

The last two expressions are simply another form of (A6-2) and A6-3). The 
following functions of u are now defined: 

Uj = cos u - K(l-cos u) (C2-35) 
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u. 


u„ 


U r 


U, 


u_ 


U r 





r (t ) 

— — — [cr(l-cos u) - sin u] 


r (t ) 
o o 


sm u 


r (t ) 
o o 

na 


[cos u - cr sin u] 


(1-cos u) sin u 

n r JU 

o o 


(1-cos u)‘ 
2 

n a 


r 2 (t ) 
o o 


t na( 1 - cos u) u. , , 4 

. 2 4 , a 1-cos u) 

in ii 4- ■ ■ 4- — ■> ■ ■ ■ 1 


sm u + 


r (t ) 
o o 


r (t ) 
o o 


- U r 


n + ( 1 - cos u) 


u 2 u 9 /n 


Uj u^/ (na) +' B RV (t Q , t) u 3 
u 2 u 5 / (na) + B RV (t o> t) u 4 
u 1 u ^/ (na) + B vv (t o , t) u 3 

u 2 u 6 /(na) + B w (t o ,t) u 4 


(C2-36) 

(C2-37) 

(C2-38) 

(C2-39) 

(C2-40) 

(C2-41 ) 

(C2-42) 

(C2-43) 

(C2-44) 

(C2-45) 

(C2-46) 

(C2-47) 

(C2-48) 
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(C2-49) 


Dj = u 1 u 7 /(na) + D RV (t Q> t) u 3 

°2 = u 2 U 7 /(na) + D RV (t o’ t) U 4 
°3 = U 1 u g /(na) + D VV (t o’ fc) U 3 
D 4 = U 2 u 8 /(na) + D VV (t o' t)u 4 


(C2-50) 

(C2-51) 

(C2-52) 

(C2-53) 


After some manipulation (C2-31) and (C2-32) become 


B (t ,t) 
r o 


D r <t o ,«) 

+ °3 V'o> io (t o )T + °4 W io (t o )T (C2 - 55) 

In these expressions for and the outer product matrices are functions 

only of t and are therefore constant for fixed t , The scalar coefficients 
o o 

are explicit functions of u (or of t through (C2-1)). Thus, as t varies only 
the scalar coefficients need to be evaluated. 


B I + B, r (t ) r (t ) + B, r (t ) v (t ) 

u 1 — o o — o o 2 — o ' o — o ' o ' 


+ B v (t ) r (t ) T + B . v (t ) v (t ) T 
2 -o' o — o o 4—0 0 —0 0 


(C2-54) 


D„ I + D, r (t } r (t ) + D, r (t ) v (t )' 

U 1 — o o — o o 2—0 0—0 0 


C3 DEFINITE INTEGRALS 

The four integral constants of the outer solution, (All-86), (All-87), (All-139), 
and (All- 140), can be written in the form 


f tk 

-10k (t k’ t o ) = J t [ B(t k’ T) ^l (T) + IlO (t k' T) ] dT 


(C3-1) 
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/ 



(C3-2) 


— llk (t k’ t o- ) 


f t k 

J 


t) 


+ Ill (l k' 




d : T 


— 20k (t k 


rt k 

•‘o 1 = J [B<t k . T >£ 2 < T) + L 2 o<‘ k . T) J 


dT 


— 21k (t k 


r fc k 

* t o ) = J [D(t k ,T ) I 2 (t) +I 2 i (t k . T) ] dT 


(C3-3) 


(C3-4) 


where 


ilO (t k- T) 


ill (t k- T) 


i20 (l k' T) 


-21 (l k’ t) 






(t k » t) - B (t k » t) b iS (t) 

(t k . T ) - D ( fc k - T ) £is ( t ^ 

(t k , t) - B(t k , t) J 2g (T) 
(t k> t) - D(t k , t) F 2S (t) 


(C3-5) 

(C3-6) 

(C3-7) 

(C3-8) 


Then's and F's are defined in. Section A1 1. 

The- integrals (C3-1) - (C3-4) cannot be integrated in closed form and must 
therefore be evaluated numerically. The integrands are finite over the. entire 
range of t giving each integral a finite value. At t = t^, however, both the 
F's and the F s are singular and the value of the integrand, which is actually 
the difference between two large numbers, must be determined by a limit 
process. This problem is eliminated by using a numerical technique such as 
Gaussian quadrature, which does not require the value of the integrand at 
either endpoint of the interval’ of integration. Thus the integrand is evaluated 
only at interior points, t o <T<t k , where both the F's and I_'s as well as the 
total integrand are finite. 
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C3. 1 Change of Independent Variable 

The"" accuracy of the integration technique for a fixed number 'of sub- intervals 

3 

within the interval t 5 t < t^ can be improved by introducing the eccentric 
anomaly through Kepler's equation 

E - e sin E (C3-9) 


(1-cos E } dE 
n 


nr = 
Then 

dr = 


r 



na 


= - — du 

na 


(C3-10) 


where (C2-1) has been introduced in the last step. Introducing (C3-10) into 
(C3-1) - (C3-4) gives 


— 10k (t k’ V 


f 


B r (t k' T) £l <T) + 


r oilO (t k’ T) 


na 


du 


(C3- 11) 


/ u 


r c — 1 1 ^k’ 

D r (^k* T > — 1 ( T ' + ‘a 


du 


(C3-12) 


— 20k (t k’ t o ) 


/ u 

° 


r n i?n (t v' T ) 

B r«k' T > H k 


du 


(C3~ 13) 


— 21k (t k 


/ u 


D r (t k’ T) —2 ^ t) + 


r o-21 (t k ,T] 

na 


du 


(C3-14) 
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whe re 


t = t, 


- j\a - e(l-cos u) sin - e sin u cos E^J 


= a £l - e cos(E^ - u)J 


(C3-15) 

(C3-16) 


u = E r E (C3-17) 

and 

u = E. - E (C3-18) 

o k o 

Equations (C3-11) - (C3-14) have an additional advantage over (C3-1) - (C3-4) 
which is not necessarily one of accuracy. It is that r appears both expli- 
citly as r^ and implicitly in determining F^ and F,,. With u (or E) as indepen- 
dent variable, r (u) can be evaluated directly. With t as independent 
o 

variable, ( T ) requires an iterative solution of Kepler's equation. Therefore 
using u as the independent variable reduces the computational requirements. 


C3. 2 Analytical Approximation for First Order Integrand 

Even though Gaussian quadrature does not require the values of the integrands 
at the endpoints, the point nearest u = 0 in (C3-11) - (C3-14) at which the 
integrands are evaluated may cause computational difficulties. This is 
because the integrands represent the differences of two large numbers when 
u is close to zero. 

til 

Theoretically the zero order outer solution passes through the center of the, 

th ^ 

k body at u = 0 (i. e. , at t = t^). Actual numerical solutions of the La'mbert 

problem may not satisfy this requirement at t = t^ but may contain a small 

residual difference between r and p. at t = t. . Experience has shown that 

-o -k k 

this small residual is greatly magnified when taking the difference between 
two large numbers. Therefore it is advantageous to have the first quadrature 
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point as far from u = 0 as possible. This may be accomplished by breaking 
the integration interval into two parts, i. e. , let each of the integrals have 
the form 


'o' 




whe re 



(C3-19) 


(C3-20) 


(C3-21) 


and 

t - t << 1 (C3-22) 

S S zC 

Over the interval from t g ^_ to t^ the integrands can be replaced by analytical 
approximations obtained from Section All. First recall that the force F^ can 
be divided into a singular part, Fj g , and a non- singular part, F^. Near t^ 
the non- singular force can be expanded in a Taylor series 


— ln^ T ' = 5:in"k> + °< T -'k> 

(C3-23) 

n, using (A3 -28) and (A3 -30) 


B(t k ,T) F ln (r) = F ln (t k ) (t k -T) + o((T - t k ) 2 ) 

(C3-24) 

D(t k' T) -ln <T> = ^ n <'k» + °' T - 'k> 

(C3-25) 
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Substituting (C3-24) and (All-45) into (C3-1) and integrating from t to t. 

S iC xC 

gives - 

/' k [B(t k ,T)F 1 (T) + I 10 (t k ,T)]dT = P 2k (t 6k - t k ) +(P4 k V 
fc sk 

+ £ln< t k>J (t 8k> t k> 2/2 

4 °[<‘sk - V 3 ' •*< t .k- t k> 2 ] (C3 - 26 > 

Substituting (C3-25) and the derivative with respect to t of (All-45) into /- 
(C3-2) and integrating gives 


r^k 

J [D(t k ,T)Fi(T) + i ll( t k .T,jdT 


'sk 


- (^k 4 ^'k 1 

4 (t sk ' V 

4 °[<‘.k - ‘k' 2 ] 


(C3-27) 


A similar analysis of the second order integrals shows that over the interval 

t , < t < t. both contribute terms of order 
S K .K 

Therefore (C3-11) - C3-14) can be written 


<l S k - V 1<>g |*8k * ‘k • 


— 10k (t k' ‘o* 


/"I 

u 


B r (t k , T ) + 


r oko ,t k' T) 


na 


du 


+ &2k (t sk " ‘k* + [&-3k + ^■4k (fc k ) + — ln^k^) 
‘ (t sk " t k ) /2 + ° | (t sk “ t k ) * ^ (t sk ~ t k ) ] 


(C3-28) 
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r U ° f r I (t, , t ) 

JiiA'V = J [ D r (t k' T) £i ,T) + 


du 


" —3k + 3 ^4k (t k ) + — ln (t k ) (fc sk “ t k ) 


+ °( (t sk - l k' 2 ) 


(C3-29) 


— 20k (t k’ t o ) 


r\ 

u • 


„ „ , „ , , . r o — 20^k’ t) 

B r (t k , T ) F 2 (t) + - 


du 


+ °(<‘sk - ‘k> lo6 l‘«k- l kl) 


(C3-30) 


f[- 


W'k-'o* = / | D r (t k’ T) L (r > + — na (t ' k ’ ~ 


du 


+ °Kk - l k> log l t »k ■ ‘kO 


(C3-31) 


The interval (t g ^ - t^) must be chosen to make the effective error of each 
approximation smaller than order p 2 . Letting 


fc sk “ fc k = ° (m ”> (C3-3Z} 

the effective error of each integral’ is 


~F 

IK 

h— < 

O 

1 

o(p 1+3n , p 2+2n ) 

^ —20 


Z+n , n\ 

a log p J 

K —a - 

o(t* 1+2n ) 

^ —2 1 

-°(> 

Z+n , n\ 

x log p j 
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2 

These error are all smaller than 0(|i ) if n >1/2. In actual numerical studies 
the value used was 

|*sk- l kl * ^ 2/3 < C3 - 33 > 

This eliminated the problems associated with the residual error in the 
Lambert solution. If only a first order solution is desired then any n >0 will 
suffice, n = 1/3 would be a likely choice. 

C3. 3 Analytical Approximation for First Order Solution 

The second order integrals are functions of the force F_ which, according to 

2 

(A5-10) is a function of r^ and r^. When integrating over the interval 

t > t > t, , r, (r) must be determined from (A6- 1 1) which itself involves an 
o k - 1 

integral function. Thus the second order integrals are actually double inte- 
grals and difficult to evaluate efficiently by numerical techniques without 
using a large number of quadrature points. 

An alternative approach is to approximate r^ with an analytical function. 

Such an approximation, valid for t near t^, is given by (All-78). The range 
of this approximation can be extended by adding additional terms and fixing 
the coefficients to give r^ the correct magnitude and slope at t = t Q . There- 
fore (All-78) is rewritten as 

£i (t) := £io (t) ' £ik log Q k (t - fc k> + £ik [ log Q k (t o - t k ) + 1 ] 

+ — 10k (t k’ V + ^2k (t " t k ) log Q k (t " t k ) + [^lk /(t o " V 

+ — llk (t k’ t o ) " ^2k [ log Q k (t o ' fc k'^ + *]] (t ~ t k ) 

- 3P 4k (t k ) (t - t k ) 2 log Q k (t - t k ) + [p 3k + 3P 4k (t k ) 

+ 6^ 4k (t k ) log Q k (t Q - t k ) + G k K 1()k (t k , t o ) 

'+ P|] (t - t k ) 2 /2 + ^ik (t - t k ) 3 + i 2k (t - fc k )4 (C3-34) 
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where 

r 10 (t) = A(t,t o ) rj(t o ) + B(t,t o ) Vj^) (C3-35) 

Taking the derivative with respect to t gives 

-l (t) = -10 {t) “ £-lk /(t “ t k ) + ^-2k [ log Q k (t ' t k ) + J ) 

+ ^lk /(t o " t k ) + — llk (t k’ t o ) " £-2k [ log Q k (fc o " + 2 J 

- 3P 4k (t - t k ) [2 log Q k (t - t k ) + ij + [i 3k + 3£ 4k (t k ) 

+ 6 £4k ( y log Q k (t o - fc k } + G k — lOk^k’ *'o ) + — k] (t " fc k> 

+ 3^ lk (t - t k ) 2 + 4^ 2k (t - t k ) 3 (C3-36) 

where 

v 1Q (t) = C(t,t o ) £l (t o ) + D(t, t Q ) (C3-37) 

At t = t (C3-35) and (C3-37) reduce to 

(C3 - 38 > 

^0 (i o> = Xlf‘o> (C3 - 39 > 


Putting t = t Q in (C3-34) and (C3-36) and solving for _£j.k anc ^ &2k S* ves 


^lk = 
^2k = 


4 ^lk jfk 

(t o - ‘k> 3 ’ <‘o ; V 2 

j|k . 3 ^ik 

“0 - l k » 3 ’ <‘o - */ 


(C3-40) 


(C3-41 ) 
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where 

ttk = - 2 £ik - JW'k' *o> - [WV V + P 2k ] <v v 

-[Psk+^k'V^kWW 

+ £^]< t o ' ‘ k ) 2 / 2 (C3-42) 

■2zk = " — llk (t k' t o ) “ [-^3k + G k — 10k ( V t o ) 

+ P*](t o -t k ) (C3-43) 

Equations (C3-34), (C3-35) and (C3-40) - (C3-43) give an analytical approxi- 
mation for r^ which has the correct behavior at both endpoints of the interval 
t <_ t t^. It is in these regions that the second order integrals are most 
susceptible to errors in r^ and therefore the effect of errors in the middle of 
the interval is minimized. 
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